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Abstract 

We explore the limits on neutrino mass which follow from a study of the 
long-range forces that arise from the exchange of massless or ultra-light neu- 
trinos. Although the 2-body neutrino-exchange force is unobservably small, 
the many-body force can generate a very large energy density in neutron stars 
and white dwarfs. We discuss the novel features of neutrino-exchange forces 
which lead to large many-body effects, and present the formalism that al- 
lows these effects to be calculated explicitly in the Standard Model. After 
considering, and excluding, several possibilities for avoiding the unphysically 
large contributions from the exchange of massless neutrinos, we develop a 
formalism to describe the exchange of massive neutrinos. It is shown that the 
stability of both neutrons stars and white dwarfs in the presence of many- 
body neutrino-exchange forces implies a lower bound, m ^ 0.4 eV/c^ on the 
mass m of any neutrino. 
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I. INTRODUCTION AND SUMMARY 



There is considerable theoretical and experimental interest at present in the question of 
whether neutrinos have mass On the theoretical side the existence of massive neutri- 

nos would require a modification of the standard SU(2) x U(l) model of the electroweak 
interaction, and could also provide some insight into the question of how other particles ac- 
quire their masses. These and other related theoretical questions, such as the missing mass 
problem in the Universe have served to stimulate continuing experimental searches for 
evidence of a nonzero neutrino mass. To date there is no compelling direct evidence for 
massive neutrinos, and the existing upper limits for Ue, I'fj,, and Ur (or for the corresponding 
antiparticles) are 

m(z/e) < 7.0 eV; m(z/^) < 0.27 MeV; m(z/^) < 31 MeV. (1.1) 

Notwithstanding the bounds in Eq. ( |1.1| ), there is mounting indirect evidence for massive 
neutrinos from a number of different experiments as discussed in Refs. [0-0. We will return 
to Eq. (0) in Sec. VIII below. 

The object of the present paper is to develop a new technique for studying neutrino mass 
which exploits the fact that massless neutrinos give rise to long-range forces. Recent interest 
in the question of weak long-range forces [p!^ -|l^ has focussed on the possible existence of 
new ultra-light bosonic fields as the mediators of such interactions. It is well known, however, 
that long-range forces can also arise from the exchange of neutrino-antineutrino {uu) pairs, 
and the history of various attempts to calculate the i/z/-exchange force is summarized in 
Ref. [|15|. The first correct calculation of the 2-body z/i/-exchange potential was carried out 
by Feinberg and Sucher (FS) [|1^ who used an effective low-energy 4-fermion interaction 
involving only charged currents. They found for the potential energy V}^^ (r) describing the 
interaction of two electrons via the i/z/-exchange diagram shown in Fig. 1, 

Kfir)=Gl/A7r\\ (1.2) 

where r = |r*i— 'r2| is the separation of the electrons, and Gp = 1.16639(2) x 10"^ GeV^'^{hc)^ 

3 



is the Fermi decay constant. Subsequently Feinberg, Sucher, and Au (FSA) [|T6[ recalculated 
Vj^\r) in the framework of the Standard Model and obtained 

V}^\r) = GU^sin'Ow + l/2)V47rV, (1.3) 

where Oyy is the weak mixing angle, with siia^ 9w = 0.2319(5). The result in Eq. ( |1.3| ) has 



been rederived recently by Hsu and Sikivie using a different formalism from that of 



FS. Numerically, V}^\r) = 3 x 10~®^r~^ eV where r is in meters, from which one can see 
that the interaction energy of two electrons arising from neutrino-exchange is smaller than 
their mutual gravitational interaction for r ^ ro = 5 x 10^'^m. However, for r ^ tq, the 
electroweak interaction arising from 7- and Z°-exchange would dominate, and hence there 
appears to be no distance scale over which V^^\r) leads to detectable effects. As we note 
in Appendix A, the neutrino-neutron and neutrino-proton coupling constants, a„ and 
respectively, are even smaller than the neutrino-electron coupling constant a^, so the same 
conclusions would hold for the corresponding 2-body potentials V^'^{r) and V^p-'(r). In what 
follows it will be shown that although the 2-body potential energy is indeed extremely small, 
the many-body interaction energy arising from neutrino-exchange can be extremely large in 
neutron stars and white dwarfs. This observation eventually leads to the conclusion that 
there is a lower bound on the mass of any neutrino or antineutrino, as given in Eq. (^.12|) 
below. 

Many-body effects arising from neutrino-exchange have been considered previously by a 
number of authors |[T8|- pO|] , often in connection with attempts to understand the gravitational 
interaction. Feynman pO[ considered the many-body interaction that would describe the 



effective 2-body force that arises when two test masses interact via neutrino-exchange in 
the presence of distant matter. Although Feynman's attempt to explain the inverse-square 
law of Newtonian gravity in terms of neutrino-exchange ultimately proved unsuccessful, he 
made an observation about such interactions which forms part of the basis of the present 
work. Feynman noted that higher order (in Gp) many-body interactions could be important 
because they depended on higher powers of the masses of the interacting objects. This is 



a special case of an observation made earlier by Primakoff and Holstein (PH) |^ in their 
classic analysis of many-body effects in atomic and nuclear systems. PH noted that in a 
nucleus containing particles the magnitude of the total A;-body interaction {k = 2,3, ■ ■ ■). 
grows as the binomial coefficient (^^^ , 

[kj ~ k\{N-k)V ^^'^^ 

which counts the number of distinct A;-body amplitudes that can be formed from N particles. 
Subsequently Chanmugam and Schweber suggested that many-body electromagnetic 



effects could be important in white dwarfs due to the presence of the binomial coefficient. 
We will return shortly to discuss the combinatorics of many-body diagrams using Eq. ( p.. 41) . 
Subsequently Hartle [jl8[ used the same phenomenological weak coupling employed by FS 



to evaluate the 0{Gp) 4-body contribution corresponding to the diagram in Fig. |^ and 
obtained the result quoted in Eq. ( |3.5CI| ) below. Using this result Hartle demonstrated 



that the energy for two electrons interacting with each other and with a spherical shell of 
matter varied with the separation r of the electrons approximately as 1/r. For the systems 
that Hartle considered the effects arising from the 4-body potential were too small to be 
of interest. Moreover, even though electrons can interact via long-range neutrino-exchange 
forces, these cannot be used to measure the electron-number of a black hole ( "a black hole 
has no neutrino hair" ) . Although these results appear to support the conclusion derived 
from the 2-body potential that neutrino-exchange forces are unimportant, we show in what 
follows that this is not necessarily the case when computing the self-energy of a compact 
object such as a neutron star or a white dwarf. In fact the self-energy of such an object 
arising from the exchange of massless neutrinos can be catastrophically large, and this result 
eventually leads to the conclusion that neutrinos cannot be massless. 

Since the large neutrino-exchange energy-density in a neutron star arises from many- 
body interactions among neutrons, it is appropriate to ask why many-body effects play 
so important a role here, but are relatively unimportant in most other circumstances. The 
explanation can be found in Table | which compares neutrino-exchange to other known forces 



with respect to three conditions which determine when many-body effects become significant. 
In order for there to be an enhancement effect arising from the binomial coefficient (^^^ , it 
is necessary in the present context to have a large number of particles interacting with 
sufficient strength in a small volume. 

The first condition is that the force be long-range, which is what allows a given particle 
to interact with many other particles in the first place. We see from Table |l| that both 
the strong interaction and the weak interaction mediated by Z^-exchange fail to meet this 
condition. In the case of strong interactions it is well known that the forces among nucleons 
saturate, i.e., that each nucleon in a nucleus interacts with only a limited number of other 



nucleons p^ , p5| . This is evidenced by the fact that the strong interaction binding energy 
of a nucleus does not increase with baryon number A as (^^^ = A{A — l)/2, as might be 
expected from a 2-body interaction, but grows approximately as A ioT A > 4. In nuclei this 
can be understood in terms of exchange forces and the existence of a repulsive hard core, 
but for a macroscopic object the finite range of nuclear forces is also significant. 

The second condition in Table | is that the force couple to a charge capable of having a 
large expectation value in a macroscopic system such as a neutron star. Evidently a long- 
range force will not produce a significant many-body interaction unless the constituents that 
are capable of interacting with one another have a (large) net value of the appropriate charge 
and a non-zero charge density. For the electromagnetic interaction a number of independent 
arguments show that the net charge Z of a typical neutron star is less than 10^^ Cq, where Cq is 
the charge of the electron, and may be as small as IO^^Cq [0. The latter estimate is smaller 
than by a factor of order 10~^°, which by itself would be sufficient to suppress many- 
body effects. However, there are in addition the suppression effects originally discussed by 



Primakoff and Holstein and these, along with the small net Z, suggest that many-body 
electromagnetic interactions should be small in neutron stars. 

The last condition in Table | requires that the interaction be sufficiently strong so that, 
when many-body interactions are in fact possible, their net effect be large. This condition 
can be understood with reference to the gravitational interaction as follows: We note that the 
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quantity which determines when gravity is sufficiently strong that many-body interactions 
would be important is the dimensionless potential $, 

where Gn is the Newtonian gravitational constant, M is the mass of the neutron star, and 
R is its radius. Using the numerical results in Sec. V below we find that $ ~ 0.2 for the 
binary pulsar PSR 1913+16, which is the "typical" neutron star we are considering. The 
fact that $ must be small is obvious since $ > 1/2 would correspond to a black hole. It 
follows that for ordinary matter many-body gravitational effects will never dominate over 



the two-body interaction, although they may lead to detectable effects in some systems 
Another way of understanding this result is to note that in natural units, where Gp and G^ 
have the same dimensions, 

Gn/Gf ^ lO-'l (1.6) 

If follows from Eq. ( |1.6D that if one were to reduce Gp by a factor of order 10'^^ so as to 
make the weak neutrino-exchange interaction have a strength comparable to that of gravity, 
then many-body neutrino-exchange effects would become relatively unimportant, just as 
gravitational many-body effects are. 

The preceding discussion can be summarized as follows: We see from Table | that each of 
the four known fundamental forces fails to satisfy at least one of the conditions that must be 
met for many-body effects to be important in a macroscopic system such as a neutron star. 
The only known interaction which is of long-range, and where there is a strong coupling to 
a charge for which a neutron star is non-neutral, is the force arising from neutrino exchange. 
It is for these reasons that many-body exchange effects can be significant for this interaction, 
while they are relatively unimportant for the others. 

To understand quantitatively how many body neutrino effects can give rise to a catas- 
trophically large energy-density, we follow the discussion of Feinberg and Sucher |]TB| who 



note that the functional form of Vjp{r) can be inferred on dimensional grounds. Since 



the Standard Model is renormalizable, the only dimensional factors upon which the static 
spin-independent potential can depend are Gp and r. Evidently the exchange of a single 
uV pair must be proportional to G'jp, from which it follows that V^^^ oc C^p/r^ in agreement 
with Eq. ( |3.19| ) below. An analogous argument shows that for A; > 3 the fc-body neutrino- 
exchange contribution W^'^^ to the binding energy of a neutron star must be proportional 
to {G^p / B?^^^){^j^ , where R is the radius of the neutron star, which is assumed for present 
purposes to contain only neutrons and to have a uniform density. Since ~ N^/k\ for 
A; < iV, it follows that 

iy(^)^lir^V. (1.7) 

For a typical neutron star (GpN/R^) = (9(10^^) [see Sec. V below], and hence it follows 
from Eq. (|1.7|) that for k <^ N higher order many-body interactions make increasingly larger 
contributions to W^''\ It can be shown that W^^^ would exceed the known mass-energy of 
a typical neutron star, and that W = Y.k W^''^ can exceed the total mass-energy of the 
Universe, as we discuss below. 

The calculation of the self-energy of a neutron star or white dwarf arising from neutrino 
exchange closely parallels the calculation of the electrostatic energy of a spherical charge 
distribution such as a nucleus arising from photon exchange. For this reason, we begin in 
Sec. II by reviewing the derivation of the familiar electrostatic (Coulomb) result, which is 
given in Eq. (|2.1| ). For present purposes the electrostatic derivation can be viewed as pro- 
ceeding in three steps: 1) First the 2-body potential is derived from covariant perturbation 
theory. 2) The 2-body potential is then integrated over a spherical volume to obtain the 
contribution from a single pair of charges. 3) The result of the preceding calculation is then 
multiplied by the factor (^'^^ = N{N — l)/2 which represents the number of pairs that can 
be formed from N charges. Particular attention is devoted to the question of integrating 
the 2-body Coulomb potential over a spherical volume, which is carried out in two different 
ways. One method uses a geometric probability technique, which is later generalized to 
apply to the many-body case. 
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Sec. Ill presents the derivation of the fc-body potential V^'^^ arising from neutrino- 
exchange, using a formahsm developed by Schwinger [p8| and originally applied to this 



problem by Hartle [18|. After reproducing the 2-body result of FS and the 4-body potential 



obtained by Hartle, the 6-body potential is derived [Eq. ( p.53|) ] and then used to generalize 



to the A;-body case [Eqs. ( p.55|) -( p.5^ )]. In Sec. IV we discuss the integration of the /c-body 
results over a spherical volume of radius R, which leads eventually to Eq. ( ^.22| ). When this 



expression is multiplied by (^^^ the result is the net /c-body contribution to the self-energy of 
the object in question. After the sum over k is carried out, the self-energy Vl^ of a spherical 
collection of neutrons can be approximated by the expression in Eq. ( |5.23| ). We then pro- 
ceed to demonstrate that, barring accidental cancellations, W/M ^ 1 for compact objects 
such as neutron stars and white dwarfs, where M is the (known) mass of each object. After 
considering in Sec. VI, and excluding, several alternative explanations for the unphysically 
large value of W/M, we turn in Sec. VII to recalculate W when the exchanged neutrinos 
have a nonzero mass m. As expected, the expression for W acquires additional factors pro- 
portional to exp(— mi?) which suppress the neutrino-exchange contribution. One can then 
calculate the minimum value of m required to reduce ly to a physically acceptable value, 
and the result is presented in Sec. VIII. From Eq. ( |8.12| ) we find for m, 



m 



> 0.4 eV/c2, (1.8) 



which is consistent with the upper limits quoted in Eq. ( p. . 1|) . 

Appendix A contains a summary of our notation and metric conventions, as well as a 
discussion of the neutrino couplings in the Standard Model. In Appendix B we present 
Hartle's (unpublished) derivation ||2^ of the Schwinger formula for W, along with some 



additional technical details of its application to the present problem. Appendix C contains a 
detailed discussion of the geometric probability formalism used in integrating the many-body 
potentials over a spherical volume. In Appendix D we discuss the formalism for summing the 
many-body contributions, and in Appendix E we generalize the Schwinger-Hartle formalism 
to apply to massive neutrinos. 



II. ELECTROSTATIC ENERGY OF A SPHERICAL CHARGE DISTRIBUTION 



As noted in the Introduction, the present calculation of the energy- density of a spherical 
neutron star arising from neutrino-exchange parallels that of the electrostatic energy of a 
spherical charge distribution. For this reason we review the derivation of the familiar result 



for the Coulomb energy Wc of a spherical nucleus [30 



^c = lz{Z-lf^, (2.1) 

where R is the effective Coulomb radius, Zcq is the nuclear charge, and el/hc ~ 1/137. 
(We use Co to denote the charge of the electron to avoid confusion with the transcendental 
number e: Ine = 1.) As noted in the Introduction, we can view the derivation of Eq. ( p.l|) 
as proceeding in three steps: (1) Determine the 2-body interaction Vci'ru) between two 
point charges; (2) calculate the interaction energy Uc between the two point charges which 
are assumed to be uniformly distributed inside a sphere of radius R; (3) generalize the 2- 
body result Uc to the Z-hodj electrostatic energy Wc by incorporating the appropriate 
combinatoric factors. 



A. Calculation of the Two-Body Potential Vc 

The electrostatic potential energy is dominated by the two-body contribution arising 
from the one-photon-exchange amplitude. In this case the evaluation of the expression for 
the potential energy Vciru) of two interacting particles is trivial, and the result is the 
familiar Coulomb potential 

Vciru) = ^, (2.2) 

where ri2 = \fi — is the distance between the charges. Although there are also in 
principle contributions from many-body electromagnetic interactions, these are small for 
reasons discussed originally by Primakoff and Holstein By contrast the weak energy 
arising from neutrino-exchange is dominated by many-body interactions, and hence the weak 
potential is correspondingly more complicated. 
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B. Integration Over a Sphere: Uc 



This is the most difficult of the three steps in the weak-interaction case. For the Coulomb 
interaction in Eq. ( p.2|) we are interested in evaluating the energy of a single pair of charges 
1 and 2 (ei = 62 = eo), having a uniform probability distribution in a volume (4/3)7ri?^, so 
that the effective number density p produced by each charge is 

From Eq. ( [2. 21 ), the potential energy arising from the Coulomb interaction between the 
charges pid?Xi and p2d?X2 centered at ri and r2 respectively is 

dUc = {pid?xi){p2d'x2)Vc{n2). (2.4) 



The average electrostatic energy Uc is then obtained by integrating Eq. (|2.4| ) over the entire 
sphere: 



= J pid^xi J p2d^X2Vc{ri2 
= elp^ dr2 rl dn rl j dVti j 



R fR r f / ^ 

■ ' ' ■ ' ' 'dn2 (— 

Vri2 



where ri = |ri| and r2 = |r*2| are measured from the center of the sphere, and dQ = dcos 9dip. 
The result in Eq. ( |2.5| ) represents the average interaction energy of a single pair of charges 
having a uniform probability distribution in a spherical volume of radius R. 

Although the 6-dimensional integral arising from the 2-body potential energy in Eq. ( p.2|) 
is straightforward, its generalization to the many-body potentials that arise from neutrino- 
exchange becomes increasingly cumbersome. We will return below to discuss an alternative 
method for evaluating the integral over a sphere, which generalizes more naturally to the 
many-body case. 
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C. Combinatorics: Wc 



As noted above, the result in Eq. ( |2.5| ) gives the energy Uc for a single pair of charges. 
For a nucleus containing Z charges the number of pairs that can be formed is Z{Z — l)/2, 
so that the final expression for the total energy Wc is 

Wc = \Z{Z - \)Uc = \Z{Z - 1)|, (2.6) 

in agreement with Eq. ( |2.1| ). For application to the k-hody problem {k = 2,3,4,...) the 
combinatoric factor Z{Z — l)/2 generalizes to the binomial coefficient ||2l| defined in Eq. ( |1.4|) 
which counts the number of combinations of k objects that can be formed from objects. 
Evidently 

which reproduces the result in Eq. ( p.6|) . 

As we noted in the Introduction, the fact that the net fc-body contribution is proportional 
to the binomial coefficient (^^^ was pointed out by Primakoff and Holstein in their classic 



paper |2T[ on many-body electromagnetic forces. Subsequently Chanmugam and Schweber 



1 2^ re-examined the question of many-body forces in electromagnetism. They observed that 
since the binomial coefficient for small k grows as 



ykj-kr (^■«' 

this factor may make many-body {k > 2) electromagnetic effects important in white dwarfs, 
where the electron density can be high. For neutrino-exchange amplitudes the presence of the 
binomial coefficient in the final expression for the interaction energy is what eventually leads 
to the large neutrino-exchange energy-density referred to in the Introduction. Specifically, 
each neutron in a neutron star carries a non-zero weak charge which couples the neutron to 
the current, so that for a neutron star A^ is of order 10^''. By way of comparison, the 
electric charge Z of a neutron star is at most of order 10^^ |^ , as noted in the Introduction. 



It follows that the enhancement of the many-body amplitudes arising from the binomial 
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coefficient is a much more significant effect for neutrino exchange than it would be for the 
corresponding electrostatic case. 

In the ensuing discussion it will be helpful to recall that the binomial coefficient (^^^ is 
a decreasing function of increasing k for fixed A^: From Eq. ( |1.4| ), 

(fc+i) _N-k 



(2.9) 



k + 1 

Hence the fc-body interaction is proportional to a combinatoric factor which, although large, 
is decreasing monotonically as k increases. From the discussion in Appendix B we see that 
the weak interaction energy arising from neutrino-exchange can be expressed as a sum over 
fc-body contributions W^^'^ as in Eq. ([Bl|), each of which will be proportional to (^^^. The 
fact that the ratio in Eq. ( |2.9| ) is less than unity when k > {N — l)/2 helps to explain why 
J2k W^'^'^ is dominated by a few terms with A; ~ A^, as we discuss in Sec. V below. 



D. Alternate Method of Integration 

Having outlined the steps that lead to Eq. (|2.1|), we return to discuss an alternate method 
for evaluating the integral of the Coulomb potential in Eq. ( p.2|) over a spherical volume. Let 
7^3 (r) denote the normalized probability density for finding two points randomly chosen in 
a uniform 3-dimensional sphere to be a distance r = ri2 apart. As discussed in Appendix C, 
the average value (g) of any function g{r) taken over a 3-dimensional spherical volume is 
then given by 

(g) = / drVsir)gir), (2.10) 







where 

r2R 







drV3ir) = l. (2.11) 



The functional form of V^lr), and its generalization Vn{f) for an ra-dimensional ball of 
radius R, are discussed in Appendix C. We henceforth drop the subscript 3 when working 
in 3-dimensions {V = V3). From Eq. ( |C7|) , with r = 2Rs, 
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Since 2R > r > 0, it is sometimes convenient to introduce the scaled dimensionless variable 
s defined above which satisfies 1 > s > 0. Plots of V{s) and its derivative V'{s) are given 
in Fig. Eq. ( p.lO|) can be viewed as an analytic Monte Carlo calculation of (g) . 



Returning to the Coulomb problem we wish to calculate (eg/r) by this method. From 
Eq. (CTD , 



el{l r) = ei dr \ — -H ^ - = — -, (2.13) 

in agreement with Eq. ( p.5|) . The advantage of this approach is that, where applicable, it 
replaces the 6-dimensional integral in Eq. (|2.5| ) by the 1-dimensional integral in Eq. ( p. 131) . 
Other useful results are given in Appendix C. 



E. Application to Massive Electrodynamics 

We conclude the previous discussion by using Eq. (|2.10|) to calculate the electrostatic 
energy of a spherical charge distribution in massive electrodynamics, i.e., for the case of 
a photon with a non-zero mass fi. This calculation completes the analogy between the 
electrostatic energy arising from photon exchange, and the weak energy arising from neutrino 
exchange when m ^ 0. When fi ^ the Coulomb potential in Eq. (|2.2| ) is replaced by the 
Yukawa potential 

VV(ri2) = ^e-^'^-, (2.14) 
and the result in Eq. (|2.13| ) generalizes to 

eg(e-^7r) = l'^F{fiR), (2.15) 

where 
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15 



1 2 1 

+ -^-^ + 



{jjRY 



(2.16) 



The "form factor" F{fiR) incorporates all the modifications arising from 7^ 0, and has the 
property that -F(O) = 1. 

We note in passing that the result in Eqs. ( 2.15 ) and ( 2.16| ) is also of interest in connec- 
tion with recent work setting limits on new forces co-existing with electromagnetism pi 



Although limits on the photon mass are quite stringent, the corresponding limits on a new 
vector force co-existing with electromagnetism are less restrictive. If eo in Eq. (|2.15| ) were 
replaced by the corresponding unit of charge / for the new field, then Eq. (|2.15|) could be 
used to set limits on and fi in appropriate systems, in a manner similar to that described 
below for the neutrino mass m. 



F. Quantum Mechanical Effects 



The preceding calculation also serves to clarify another issue which arises in neutrino- 
exchange, namely, why the charge distribution can be treated classically even in an object 
such as a nucleus or a neutron star where quantum effects are important. In particular one 
may ask what role the Pauli exclusion principle plays for the protons in a nucleus or the 
neutrons in a neutron star in the presence of a long-range force. For present purposes we 
can invoke an argument due to Fermi to demonstrate that the self-energy of a nucleus 
or neutron star arising from long-range forces can in fact be approximated by the classical 
result, as we have done. (We will return in Sec. V below to discuss the effects of the Pauli 
principle for the neutrinos.) 



III. THE MANY-BODY POTENTIALS ARISING FROM 
NEUTRINO-EXCHANGE 



Following the discussion in Sec. II, we present in this section a derivation of the many- 
body spin-independent potentials arising from neutrino-exchange. We will begin by using 
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the Hartle formalism [0,^ to re-derive the original 2-body result of Feinberg and Sucher 
(FS) arising from the diagram in Fig. |1| [p!5|-p!7|. It will then be shown that fc-body potentials 
where k is odd make no contribution to the neutrino-exchange energy of a spherical neutron 



star. After reproducing the 4-body result of Hartle [|I^, we derive the 6-body potential 
which is then used to infer the relevant combinatoric factors in the general fc-body potential. 
For the sake of definiteness we assume that the external fermions are neutrons, as would be 
the case in a neutron star, and we therefore suppress the corresponding subscripts on V and 
W. 

Our starting point is the Schwinger formula in Eq. (|B34|) , 

W^ = ^Tr|/_^rfE ln[l + ^iV,7,(l + 75)4°H^)]}, (3.1) 

where Tr denotes the generalized trace defined by Eq. (|B13|) , and ln[l + ■ ■ ■] represents the 
infinite series 

CO A A; 

ln(l + A) = -E(-l)V- (3.2) 

k=i ^ 

The factor 1/k can be understood as the product of l/k\ arising from the perturbation 
expansion of exp(z£/(x)) where £/(x) is given in Eq. (|B16| ), and {k — 1)! which represents 
the number of ways that the k external currents can be attached to the closed neutrino 
loop. There is an additional factor of k\ which counts the number of ways that the momenta 
carried by the external currents can be assigned to the k vertices. The net result is that the 
expansion of ln[l + ■ ■ ■] in the Schwinger formula gives rise to (A; — 1)! topologically distinct 
diagrams, but there is no overall /c-dependent numerical coefficient |3^ . For = 4 there are 
6 independent diagrams: These are given by the 3 shown in Fig. ^ along with an additional 
3 diagrams obtained from those shown by reversing the direction of the internal neutrino 
momentum, as in Fig. ^. For /c = 2, (A; — 1)! = 1 which agrees with our expectation that 
there be only one diagram in order G|., since both senses of the neutrino momentum are 
topologically equivalent. 
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A. The 2-Body Potential 



Combining Eqs. ( ^.1| ) and ( |3.2| ) the k = 2 contribution to W is given by 

xtr {7^(1 + 75)4°^ (ri2, ^)7.(1 + l5)S?{f2u E)} N^{xi)N,{x2), (3.3) 

where ru = (xi —X2), N^{xi) and A^j,(a;2) denote the external neutron currents at Xi and X2, 
and Sp\fij,E) is given by Eq. ( P42|) . The overall minus sign arises from the expansion in 
Eq. ( p.2|) , and tr denotes the trace over the Dirac matrices in {■■■}• The factors of (1 + 75) 
can be anticommutated past Sp^ and, using the relation 

(1 + 75)^ = 2(1 + 75), (3.4) 

we have 

tr {■■■} = 2tr {sPir2i, E)^^sP{fu, E)^,{1 + 75)} . (3.5) 

Since we are interested in computing the self energy of a static collection of neutrons (which 
are assumed to have no net polarization), the neutron currents are given by 

N^{xi) = ipi6^A] N^{x2) = ip2Su4, (3.6) 

where pi = p2 = p is the number density of neutrons in the neutron star, and the factors of i 
arise from Eq. (A4). Combining Eqs. (|3.5| ), (|3.6| ), and (A8) we see that the term containing 
75, whose trace is proportional to e^i^xp, makes no contribution. The remaining terms give 

X tr [7^747/374] r/,(21)A^(f2i, E)r/^(12)Ai.(fi2, E). (3.7) 

Using 

tr[7Q 747/374] = 4(25^,4^/34 - Safs), (3.8) 
the expression in curly brackets in Eq. (|3.5|) can be written as 
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{(3.5)} = 4[2r74(21)r74(12) - r/(21) ■ 7^(12)] 

= A[E^-du-d2i], (3.9) 



where di2 = d/dri2- Combining Eqs. (|3^) and ( p.7|) we have 

X {E^Apir^i, ^)A^(fi2, ^) - ■ 4Ai.(f2i, ^)A^(fi2, E)} 
— I y pid xi J p2" y 



TT 



X 



) \e'- du-d,, . (3.10) 



471 J I ri2r2i V ^12'"21 



We note from Eq. ( [B42D that the operators du and 821 act on the respective coordinates 
ri2 and r2i as if these were independent, notwithstanding the fact that ri2 + r2i = 0. This 
apphes as well to all the derivative terms that appear in the fc-body amplitudes. 

Following Hartle [29| the integral over E can be evaluated by considering the functions 
In{z) defined by 

/oo 
dE E-e*l^l("+") 
-00 

2 fn°° dE e^l^l(^+^^) = Uz) even n 

^ ' (3.11) 

odd n, 

where z = ryi + '"2i- Since \E\ is an even function of E^ In{z) is nonzero only for even values 
of n. An elementary integration gives 

Io{z) = (3.12) 



and differentiating Eqs. ( p.ll| ) and ( |3.12| ) with respect to z leads to 



Continuing in this way we find 



Combining Eqs. (|3.10|) and ( p.l4| ) allows ly(^) to be written as 



\72 



X 



1 , ^ 1 
-f \ + -Ci2 ■ C21 -, ■ ^ 

ri2r2i[ri2 + r2iY 2 ri2r2i(ri2 + r2i , 



(3.15) 



We note from Eq. ( |2.5| ) that the integrand in Eq. (|3.15| ) has the same form as in the analogous 
electromagnetic case, with the quantity in curly brackets in Eq. (|3.15|) corresponding to the 
2-body potential V'^'^\ri2). Thus 



V2 



ldi2 ■ 821 ^ 

2 ri2r2i(ri2 + rai, 



(3.16) 



7l2'"2l(n2 + ^"21)^ 

Since r = ri2 = T21, first term in square brackets in Eq. ( 3.16 ) reduces to l/8r^. In the 
second term we note that the gradients act on a function which depends only on ri2 and 
r2i, and hence we can write 



di2 ■ 821 



d 



d 



d d 



d d 



ri2- 



r2i- 



ri2 ■ ^21- 



' dru J \ dr2i J dr^ dr2i dr^ dr2i ' 

Using Eq. ( gTzD the expression in square brackets in Eq. (|3.16| ) can be written as 



(3.17) 



1 5 _ -1 



(3.18) 



and hence, 



Gpa'l 1 



(3.19) 



4^3 ^5 

Eq. ( |3.19| ) gives the original FS result |1T5|-[T7[1 when we set a„ = 1, which is the value 
appropriate to the charged-current model of the weak interaction assumed by FS. 

For later purposes it is interesting to note that the functional form of V^'^^r) can be 
inferred on dimensional grounds, as noted originally by Feinberg and Sucher |T^. The 
only dimensional quantities upon which a static neutrino-exchange potential can depend are 
Gp, r and (possibly) the masses of the external particles. However, in the non-relativistic 
limit appropriate to a static potential, bilinear covariants such as u{p')'yx{l + '~i^)u{p) are 
independent of the mass of the fermion characterized by the spinor u{p). Thus the only 
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relevant dimensional parameters are Gp and r and, since the 2-body operator is proportional 
to G^., it follows that V"*^^^(r) oc G\/r^. Implicit in this argument is the assumption that no 
other dimensional parameters are present and, since the standard model is renormalizable, 
this will indeed be the case. This argument holds even in the framework of the (non- 
renormalizable) charged- current model originally assumed by FS, since the regularization 
procedure employed by FS to extract the long-distance behavior of V^'^\r) introduces no 
additional mass parameters. 

We conclude the discussion of the 2-body potential by demonstrating quantitatively that 
its effects are too small to be detected at present in any known system. Let us consider 
the analog of Eq. ( p.l9| ) for electrons, for which the "background" gravitational interaction 
would be smallest. Reinstating h and c, and substituting = (2sin2 6'vi/ + l/2), where 

i 

sin^^VK = 0.2319(5), (3.20) 

we find 



Using 



G 

^ - 1.16639(2) X lO^^GeV"^ (3.22) 



(he) 



leads to 

V-(^)(r) = 3xl0-^ % 'f^ (3.23) 
^ ' (r/1 mf ^ ' 

The magnitude of the corresponding force Fi2{r) = —^V^'^\r) is then given by 

2 X 10"^°° 

I^i2(r)| = Newtons. (3.24) 

(r/1 mj° 

To appreciate how weak this force is, we compare |Fi2| to the gravitational force Fi2°'^ at a 
nominal separation r = 1 km, 
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iFiofr = 1 km) I 

' ^ 4 X 10"^^ (3.25) 



\Ffrir = l km) I 

As r decreases |Fi2(r)| increases more rapidly than does |-Fi2^"''|, and these forces become 
equal at r = 5 x lO^^m. However, at this separation electromagnetic and weak forces 
arising from Z'^-exchange would be much larger than I-F12I, and hence there appears to be 
no distance scale over which the presence of the 2-body neutrino-exchange interaction can 
be directly detected [0,53]. 



B. The fc-Body Contribution when k is Odd 

In this subsection we show that there is no 3-body static potential arising from neutrino- 
exchange. This result is then generalized by demonstrating that for odd values of A; > 5, 
the fc-body potential exists but averages to zero when integrated over a spherical volume. 
The net result is that the neutrino-exchange energy of a spherical nucleus or neutron star is 
given by a sum of k-hodj contributions where k is even. 

As noted previously, there are two independent diagrams for k = 3 and these are shown 
in Fig. Expanding Eq. ( p.l|) to 0{G^) we find for the contribution from diagram ^(a) 

27r y a/2 / J J-00 

X tr{7A(l + 75)4°^(12)7m(1 +75)4°^(23)7.(1 + 75)4°^(31)} 

X N^{x^)N^{x2)K{x^)], (3.26) 

where Sp\l2) = Sp\fi2, E) etc. Using Eq. ( |3.4| ) the expression in {■ ■ ■} for diagram |^(a) 
can be simplified to 

tr{diagram 5(a)} = 4 tT{-f^sP{12)-f^sP{23)-f,sP{31){l - 75)}. (3.27) 



The contribution W^, from diagram ^(b) is given by the same expression as in Eq. (p.26|) 
except that tr{- • ■} is replaced by 

tr{diagram 5(6)} = 4 tii^^sP {13)j,sP {32)j,sP {21){1 - 75)}. (3.28) 
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The familiar arguments used to derive Furry's theorem can now be adapted to 
show that the 3-body contribution vanishes. We introduce the charge-conjugation matrix C 
defined by 



C"i(l-75)C=(l-75)' 



(3.29a) 
(3.29b) 



where the superscript T denotes the transpose matrix. In the Dirac-Pauh conventions C 
is given (up to an overall phase) by C = 7274. By virtue of Eqs. ( p.29a|) and (|3.29b|) , the 
neutrino propagator in Eq. (|B42|) satisfies 



C-'sf{n,, E)C = -7'^' ■ r/(zj)A^(r-;,, E) = S'^ (-r.,, -E) ^ Hj 



(0)T. 



(3.30) 



Using Eqs. (|3.29a|) , (|3.29b|) , and ( p.30|) , Eq. (|3.28|) can be rewritten in the form 



4 tr{7A4°^(13)7.4°^(32)7^4°H21)(l - 75)} 

= 4tr{(-7r)4°^^(-13)(-7n4°^^(-32)(-7j)4°^^(-21)(l - 75)^} 

= -4 tr{7A4°^(12)7M4°^ (23)7.4°^ (31)(1 + 75)}. (3-31) 

The last step of Eq. ( ^.31| ) follows by noting from Eq. ( p.30| ) that we can replace —fij 
by fji in Sp\—ij), so that Sp\—ij) is the same as Sp\ji) except for the sign of E. 
However, by virtue of the symmetric limits of the integration over E, only even powers of 
E contribute when the products of the neutrino propagators are expanded. Since these give 
the same contributions for ±E, we can effectively replace Sp\—ij) by Sp\ji) in Eq. ( p. 311) . 
Combining Eqs. (|3.27|) and ( p.31|) we see that the terms independent of 75 cancel exactly, 
which is Furry's theorem, while the terms proportional to 75 add yielding 



|y(3) _ 1^(3) ^ ^r 



(3) 



-i ( Gpan 
2^ 



X 



//•oo 
d^xid^X2d^X3 / dE 
V — / J —00 

ti{jxSP{12)j,sP{23)j,sP{31)j,}Nx{x,)N^{x2)N,{xs) 



(3.32) 
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Using Eq.(A8) we see that the trace in Eq. (|3.32|) is proportional to the Levi-Civita tensor 



^^lu\p^ which reduces in the non-relativistic hmit to the permutation symbol eiji. It then 
follows that the 3-body potential ^^'•'^■'(^12, ris, ^23) will contain terms of the form 

(^i 12) i{rzi)j{r2z) I = ^2 ■ (^31 x ^23) = 0, (3.33) 

since fi2 + rsi + r23 = 0. Thus there is no 3-body potential arising from neutrino exchange. 

The preceding arguments can be generalized to show that one can ignore the contribu- 
tions from fc-body potentials when A; > 5 is odd. As we discuss in the following subsection, 
the k-hodj potential arising from the expansion of Eq. ( |3.1|) receives contributions from 
^(A; — 1)! topologically distinct pairs of diagrams, each being sum of two diagrams, repre- 
senting the two senses of the neutrino-loop momentum. For the j-th such diagram the sum 
of the contributions from the two senses of the neutrino-loop momentum gives rise to an 
analog of Eq. (|3.32|) which, if k is odd, has the form 



00 



27r V 72 , 

x2hT{jxSP{12)j,sP{23) ■ ■■SP{kl)j5}Nx{xi) ■ --N.ixk). (3.34) 

Evaluation of tr{- ■ •} leads again to an expression which is proportional to e^^^Ap? and hence 
to eiji in the nonrelativistic limit. For A; > 5 an expression such as Eq. ( 3.33|) need not vanish. 



since eiji can be contracted with 3 linearly independent vectors, and hence the corresponding 
potential is in general nonzero. However, the integral of any such potential over a spherical 
volume is zero: In analogy to the electrostatic case discussed in Sec. II, the integrated 
energy can depend only on G^? and the radius R of the sphere. Since a nonzero scalar 
product cannot be formed utilizing only e^i, R, and R = RR, the /c-body contribution to 
the energy of a spherical charge distribution must vanish when k is odd. This conclusion does 
not necessarily hold for an asymmetric charge distribution, but in the present context any 
contributions from k-odd terms would be proportional to the (presumably small) deviations 
of the matter distribution in a white dwarf or neutron star from spherical symmetry. 
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C. The 4-Body Potential 



Following the discussion in Appendix B we note that the expansion of Eq. ( |B34|) in 
powers of Gp leads in order G^p to the one-loop diagrams shown in Fig. ^, each having 
k vertices. Since there are {k — l)!/2 distinct ways that the integers l,2,3,...,fc can be 
arranged on the perimeter of a circle, there are altogether {k — 1)! topologically distinct 
diagrams in order k, where the additional factor of 2 takes account of the two senses of 
the neutrino-loop momentum. In each order we will refer to the diagram in which the 
vertices are labeled sequentially 1, 2, 3, . . . , /c as the "standard" diagram. When computing 
the fc-body potential, each of the (A; — l)!/2 possible sequences of integers leads to a distinct 
potential arising from the sum of the two diagrams with opposite senses of the neutrino-loop 
momentum. Consider, for example, the three topologically distinct diagrams corresponding 
to A; = 4 shown in Fig. |[ Since the dependence of any term in the potential on rij arises 
from the neutrino propagator connecting i and j, the standard diagram. Fig. 0(a), will be 
a function of only the variables ^12, "^23, ''^34, and r4i, but not of r24 or ri3, which are not 
connected by neutrino propagators. Thus the 3 pairs of diagrams represented in Fig. ^ 
give rise to 3 distinct contributions to the total potential in order G^, as we show explicitly 
below. However, on symmetry grounds each of these 3 contributions leads to the same result 
when integrated over a spherical volume. Hence in practice it suffices to evaluate the pair 
of diagrams containing the standard sequence of vertices 1, 2, 3, . . . , fc, in the anticipation 
that the integrated result will eventually be multiplied by {k — l)!/2 to obtain the total 
contribution in order k. Since any specific set of k particles can be chosen in ways from 
among particles, there is an additional factor of present in the final result, and it is 
this factor which is ultimately responsible for the large neutrino-exchange energy- density. 

We can summarize the preceding discussion and Appendix B as follows: The expansion 
of the Schwinger formula in Eq. ( P34|) leads to a set of irreducible 1-loop Feynman diagrams 
which describe the 0{G\) contribution to the fc-body potential. The /c-body potential can 
also receive contributions from diagrams containing higher powers of Gp, but these are 
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suppressed both by and by various mass factors. Hence for practical purposes we need 
consider only the one- loop diagrams arising from the expansion of Eq. ( |B34| ) . Each of these 
contains k vertices, which represent in configuration space the coordinates fi{i = 1,2, k) 
of the k particles. The k vertices are connected by k legs representing the variables Vij = 
Ifijl = \fi — fj\. There are k{k — l)/2 such variables that arise, and {k — l)!/2 topologically 
distinct diagrams, although each diagram depends on a subset containing only k of the 
variables. Notwithstanding the fact that the variables r^j satisfy a constraint of the form 

ri2 + r23 + r34-\ h ni = 0, (3.35) 

the evaluation of the fc-body potential for k > 2 should be carried out treating the fij as if 
they were in fact independent. 

We turn next to the detailed form of the 4-body contribution. It follows from the pre- 
ceding discussion that we can confine our attention to the standard diagrams. Figs, ^(a) and 
^(a'), which incorporate both senses of the neutrino loop momentum. Expanding Eq. ( [B34| ) 
to 0{Gp) and using Eq. ( pj.4D we find for the contribution from diagram Q^a), 

2tt \ v2 / "'-oo 

x23tr{7^4°)(14)7.4°)(43)7A4°n32)7.4°^(21)(l-75)}T^.A.(xi,a;2,a;3,X4), (3.36) 

where 

T^uXa{Xi,X2, X3, X4) = Nf,{xi)N^{x2)Nx{x3)N„{x4,) . (3.37) 
By use of Eqs. ( 3.29a| ) and ( |3.29b|) , tr{- ■ ■} in Eq. ( 3.36| ) can be written in the form 



tr{diagram 4(a)} = tr{5?)(12)7,4°)(23)7A4°^(34)7.4°) (41)7^(1 - 75)}. (3.38) 

The contribution from diagram §(a'), corresponding to the standard diagram with the loop 
momentum reversed, has the same form as Eq. ( p.36| ) except that tr{- ■ •} is replaced by 

tr{diagram 4(a')} = tr{4°) (12)7,4°) (23)7a4°^ (34)7^4°) (41)7^(1 + Ts)}- (3.39) 
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We see that when the contributions in Eqs. ( p.38| ) and (|3.39|) are added the terms pro- 
portional to 75 now cancel, whereas they added in the 3-body case. Similarly, the terms 
independent of 75, which cancelled previously to yield Furry's theorem, now add to give 



—i ( Gpar, 



(4) 



4 

I 2^ J (fxi...Sxij dE 



27r V V2 

xtr{Sf (12)7,4°) (23)7a4°^(34)7.4°)(41)7^}T,,,a.(xi, • • • ,^4). (3.40) 
As in the 2-body case we are interested in the static potential which arises from 

T^uXa{.Xi, . . . , X4) = (z)VlP2P3P45/x45j/45A45f74, (3.41) 

and hence the trace in Eq. (|3.4CI|) is proportional to 

tr{(3.40)} ~ tr[7 ■ 7/(12)747 ■ 7/(23)747 ■ r?(34)747 ■ 77(41)74], (3.42) 

where our notation is the same as in Eq. ( p.?]) . It is convenient to simplify the trace by 
writing 

747 ■ 7/(12) = -(7 ■ 9i2 + 74^)74 = -7 ■ ^(12)74, (3.43) 
so that tr[- ■ ■] in Eq. ( p.42|) assumes the form 



tr[- ■ ■] = tr[7 ■ 7/(12)7 ■ ^(23)7 ■ 7/(34)7 ■ 7/(41)] 

= A{E'^ - ■ 823 + 823 ■ ^4 + ^4 ■ ftti + di2 ■ 834 + di2 ■ dii + 823 ■ d^) 

+ [012 ■ d23)03A ■ ail) - (^12 ■ 83^)023 " ^41) + 0u " d^l)023 " 4)]}- (3.44) 



Combining Eqs. ( p.40| ), ( p.41| ), and ( p. 441) , we can write W^^^ in the form 



^(4) _ zi f^p^^ 



2^ J pid^xi j p2d^X2 J P3d^X3 j p4d^X3 J dE 



27r V V2 

xtr[- • ■]Ap(fi2, E)Ap{r23, ^)A^(f34, ^)A^(7='4i, E), (3.45) 



where tr[- ■ ■] denotes the expression on the right-hand-side of Eq. ( [3.441 ). Using Eq. ( [B42[ ) 
the product of the functions Apivij, E) can be written as 
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Ai.(fi2, E)Ap{f23, E)Ap{f3,, E)Ap{f,u E)=( — ) . (3.46) 

If we denote the sum of the rij and their product by 5*4 and P4 respectively, 

5'4 = ri2 + r23 + r34 + r4i, (3.47a) 

Pi = ri2r23r34rii, (3.47b) 



then the integral over E can then be expressed in terms of 5*4 by using Eqs. (p.ll|) and 
( p.l4| ), and then making the replacements 



E^ ^ 2^-^, (3.48a) 

E^ ^ -2i|i, (3.48b) 
04 

E° 2i^. (3.48c) 
04 



Combining Eqs. (|3.45|) - (|3.47b| ) we can extract the contribution to the 4-body potential V^'^^ 



from diagrams 4(a) and 4(a') by writing 

iy(^)[4(a) + 4(a')] = J pi^xi J p2d'x2 J p^d'x^ J P4d'x,V^^\fu,f23,r34,r4i), (3.49) 



where 



V'('')(ri2,f23,r34,r4i) = ( — ^ ) -^l- 



4! 



\ V2 J ^TcHP^S! 

+ 2! (J)i2 ■ d23 + 823 ■ ^34 + (934 ■ 4l + <9l2 " ^34 + 812 " fttl + ^23 " <94l) 



+ 



{di2 ■ <923)(534 ■ <94l) - (<9l2 " 534)(<923 " ^41) + (<9i2 " 94l)(923 " <934) 



-P4>S'4 
1 

P4S4 . 

(3.50) 



The expression in Eq. ( |3.50| ) reproduces the 4-body result of the Hartle |]T8[ up to some minor 



misprints in that reference. We note again that the gradient operators dij act on P4 and 5*4 as 
if all the coordinates fij are independent, notwithstanding the fact that ri2+r23+r34+r4i = 0. 



27 



D. The 6-Body Potential 



The general discussion that introduced the 4-body potential in the preceding subsection 
can be taken over immediately for the 6-body case. Our purpose in deriving the 6-body 
potential V^^^ is to explicitly exhibit the various phases and combinatoric factors in a way 
that will allow us in the next subsection to obtain the general /c-body result V^''\ From the 
preceding discussion we note that for A; = 6 there are 60 pairs of topologically distinct dia- 
grams, with each pair representing the two possible senses of the neutrino loop momentum. 
The 6 legs of the standard diagram correspond to the variables ri2, r23, r34, r45, rse, and rgi 
which satisfy 



ri2 + + + f^^ + fge + rgi = 0, 



(3.51) 



and 5*4, P4 in Eqs. (|3.47a|) and ( |3.47b|) are now replaced by 



= ri2 + r23 + r34 + r^r, + rgg + rgi, 



(3.52a) 



Pe = ri2r23r-34r45r56r6i. 



(3.52b) 



Using the formalism of the previous subsection the expression for V^^\fi2, . . . ,rei) can be 



written as follows: 




(3.53) 



where 




'j- - {ifAl di2 ■ 823 + ^34 ■ ftis + ^56 ■ ^61 + di2 ■ du + d23 ■ <9z 



'45 



+ di2 ■ ^45 + ^23 ■ ^34 " ^12 " ^6 " ^^23 " 4>1 - ^12 " 4l 



- ^23 ■ ^56 + ^34 ■ del + ^34 " d^G + 9^5 ' Qqi + ^45 " ^6 



+ (0^2! 012 ■ 923)(4 ■ 556) + (^12 ■ d23){du " del) + (4 ■ dei){d,5 ■ 4) 

+ (^12 ■ 923)034 ■ Sse) + (^12 ■ d23)0A5 " ^il) " (<934 " 9^^)045 ■ Oqi) 
- 012 ■ 934)023 ■ Oqi) + 012 ■ 034)045 ■ ^g) + (<923 " 4l)(<945 " ^e) 
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del) 
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del) 
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■ del) 
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del) 
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del] 
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■ del) 
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■ del 
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[di2 


056) 


(^23 


del) 


- (<934 


045) 


{012 


d^e) 


- (^34 


d45 


+ 


012 


945) 


(^23 


934) 


+ (^12 


945) 


05e 


del) 


+ (4 


d34] 
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023 


•945) 


(^6 


del) 


-(^12 


934) 


023 


^45) 


+ (^12 


d34] 


+ 


012 


^^23) 


(^4 


945) 


+ (^12 


823) 


05e 


del) 


+ (4 


d45] 



+ 
+ 



+ 
+ 



- {di2 ■ 

dl2 ■ <934 
dl2 ■ <945 
^23 ■ ^34 
^23 ■ <945 
^34 ■ ^45 
dl2 ■ <945 
dl2 ■ d23 



23) 
d23 ■ "961 

•923 ■ del 
di2 ■ del 
di2 ■ del 
di2 ■ del 

d23 ■ C^34 
^34 ■ ^45 



){d34 ■ dei){d45 ■ di 



'56 

(^45 ■ <956) - 
(<934 -46) + 
(^15 -46) + 

— * — * 

(^34 ■ <956) - 
(^23 ■ ^56) - 
(<956 -961) + 



<9l2 ■ f^34 

<9i2 ■ AtS 

f^23 ■ C^34 
d23 ■ C^45 
^34 ■ ^45 

<9l2 ■ f^34 



+ (<9i2 ■ d23){du ■ di 



— * 

)[d45 


del) 


)(534 


d5e) 


)034 


■ del) 


)045 


d5e) 


)045 


del) 


)012 


del) 


)034 


del) 


){d23 


■ d^e) 


){d23 




J 1^56 


^61 J 


){d5e 


del) 


)05e 


del) 


)045 


del) 



'56 

(•923 ■ d^e)04b ■ del) 
023 ■ dr>e)034 ■ del) 
012 ■ d5e)045 ■ del) 
012 ■ d5e)034 ■ del) 
012 ■ d5e)023 ■ del) 
023 ■ d45)05e ■ del) 



PeSe 



PeSe 



(3.54) 



){d5e ■ del) 

It is helpful to understand the sources of the various numerical factors in V'^^^ as a prelude to 
discussing the fc-body potential V^''\ The overall coefficient (—2/277) arises from Eq. ( |B34D , 
with an additional minus sign from the expansion of ln(l + A) in Eq. ( |3.2| ). Each of the 
neutrino propagators in Eq. ( [B4U| ) contributes a factor (i/Aii), and In{z) in Eq. ( |3.14D is 



responsible for the factor (i)^. The 6 factors of (1 + 75) eventually collapse to a single factor 
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by use of Eq. ( p.4|) , leaving a numerical coefficient of 2^. An additional factor of 2 arises 
from the diagram with the reversed loop momentum and this, along with the factor of 2 
from In{z) in Eq. ( |3.14| ), explains the factor 2^. The Dirac trace leads to the coefficient 4 
multiplying and in 6th order the overall strength of the 6-body interaction is given by 
iGFajV2f. 



E. The fc-Body Potential 



The preceding discussion can be generalized in a straightforward manner to the fc-body 
potential V^'^'^. As we discuss in Sec. V below, it is sufficient for present purposes to consider 
only the single term in V^''^ which is independent of derivatives, and hence we wish to 
exhibit this term in complete detail. We will also exhibit the coefficients of the derivative 
terms, whose forms can be obtained by generalizing Eq. ( |3.53|) in an obvious way. Collecting 
together the various powers of i and factors of 2, we can write the contribution to V^''^ from 
the standard diagram in 0{Gp) as follows: 

1 



Att) [ y/2 



where 



{i)''k\ 

+ {i)''-\k-2)\ 



Im ) 



+ ■■■ + 



Yl'^''\dab ■ drs){dcd ■ di 



Irn I 



PkSk 



(3.55) 



(3.56) 



Sk = ri2 + + H h r^i, 

Pk = ri2r23r34 ■ --rki. 



(3.57) 
(3.58) 



In Eq. ( ^.55| ) the notation • • • , Z)*''^'' symbolically represents the terms containing 

products of 2, 4, . . . , A; derivatives, which act to the right on the indicated functions of Pk and 
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Sk- The numerical coefficient of ^^'^^ can be understood in light of the preceding discussion 
as follows: The factor of 1/2tt arises from the product of the overall coefficient (—2/277) of 
W in Eq. (P34| ) , Eq. (|3.2| ), and a factor of i extracted from In{z) in Eq. ( p.l4|) , where n —>■ k, 
k — 2, k — 4:,... as appropriate in ^C^). The factor of 2'^"'"^ is the product of 2''~^ arising from 
k factors of (1 + 75), and factors of 2 from /„(-2) and from the diagram with the reversed 
neutrino-loop momentum. The overall strength of the A;-body contribution is determined 
by [G Fdnl V^Y ^ and the neutrino propagators contribute the factor [I/A-kY and the phase 
(z)'^. This phase, when combined with the k powers of i resulting from the generalization 
of T^i/Ao-(a;i, . . . , 3:4) in Eq. ( p.41| ), leads to the overall phase {iY^ = +1 for even k. As 
before, the Dirac trace contributes the coefficient 4 of and /„(-2) gives the coefficients 
k\, (k — 2)\, . . . of the individual terms in Collecting together the various factors in 

Eq. ( p. 551) we have (for k even) 

V^'\fu,...,ni) = ^ {^)''^^"^(S,,P,). (3.59) 

It is worth noting that the coefficient of ^^''^ is real, as is each of the terms in However, 
because the individual phase factors {iY, (i)^~^, . . . can (depending on k) assume the values 
(±1), various partial cancellations take place. Since the individual terms are unphysically 
large, these cancellations do not eliminate the problem arising from the energy- density due 
to neutrino exchange, but produce other novel unphysical effects, as we discuss in Sec. V 
below. We will thus be led eventually to the conclusion that because all these effects are 
consequences of the assumption that neutrinos are massless, there must in fact be a lower 
bound on the mass of any neutrino. 



IV. INTEGRATION OF THE NEUTRINO-EXCHANGE POTENTIALS OVER A 

SPHERICAL VOLUME 

Having derived the analytic forms of the fc-body potentials in the previous section, we 
turn to the problem of integrating these over a spherical volume, in analogy to the electro- 
static case. 
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A. The 2-Body Potential 



From Eq. ( ^.211 ) the 2-body potential between neutrons arising from neutrino exchange 

is 

^"'('■' - - ^- 

where r is the distance between the neutrons. In the absence of a cutoff hmiting how 
small r can become, the integral of V^'^\r) over a spherical volume would be singular. This 
contrasts with the electrostatic case discussed in Sec. II, where the radial factors in the 
volume element d^Xid^X2 in Eq. ( p.5|) offset the contribution from the potential eg/ri2 and 
lead to a well-behaved result. In practice a natural cutoff on r exists in the systems of 
interest to us, specifically in nuclei and in neutron stars. In the former case it is well known 
pO| that the nucleon-nucleon interaction has a strong repulsive component (the "hard core" ) 
which prevents the nucleon-nucleon separation from becoming smaller than approximately 
0.5 fm. The dynamics of neutron stars, although less well known, would also require a hard 
core. Combining Eq. ([4.1|) with Eqs. ( |2.10| ) and ( p.l2| ) we find for the 2-body contribution 
arising from neutrino exchange 

Jtc 

r2R 



/■2« / K 



K, 



2R^r^ \ 4i?2 J 4Rir^ y 2RJ IQR^ 



3k 1 



(4.2) 



2 i?V2' 

Eq. ( [4.2| ), which is the analog of Eq. ( |2.5| ) in the electrostatic case, gives the average inter- 
action energy for a single pair of neutrons having a uniform probability distribution in a 
spherical volume of radius R. 

It is worth commenting on the implication of the fact that the integration of V^'^\r) 
over the spherical volume leads to the introduction of an additional dimensional parameter, 
namely r^. We note to start with that since < R the presence of Vc strengthens the 
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dimensional arguments given in the Introduction concerning the approximate magnitude of 
the contributions from neutrino exchange. Secondly, we anticipate the ensuing discussion 
of the fc-body contributions by noting that for k > A the integrations of the potentials T^^^-* 
lead to expressions for U^'^'^ which are well-behaved as Tc ^ 0. Hence with the exception of 
f/^^'', which makes a negligible contribution to the total energy density, all t/^*"'-' for k > A 
depend only on Gp and -R, as noted in the Introduction. 



B. The 4-Body Potential 



The new feature which arises in the integration of V^^'^ (ri2, . . . , r4i) , as given in Eqs. ( p.49[ ) 



and ( p.50|) , is the presence of angular factors which result from the scalar products of the 



derivative terms. Although one might be tempted to assume that any such factors average to 
zero when integrated over a sphere, this is not the case as we discuss below. The evaluation of 
the terms involving derivatives can be simplified by generalizing the result in Eq. ( |3.17| ). To 



start with we recall from Eq. (|B42| ) that all the derivatives have their origin in the neutrino 
propagator S^p\ where dij acts on rij appearing in Ap{rij^E). In principle the presence of 
all derivatives can be eliminated at the outset by evaluating dijAF{rij,E) before carrying 
out / dE. Although it is more convenient to evaluate / dE first, one must recall that the dij 
which then appear in are to be understood as acting on the corresponding T>i j clS if all 
the fij were independent. It is straightforward to show, for example, that the same 2-body 
result is obtained if one evaluates dijAp first, rather than leaving the derivatives to the end 
as was done in arriving at Eq. ( p.l6| ). 

Since a derivative such as di2 in V^'^^ acts on a function /(ri2, . . .) which depends only 
on the magnitude of fi2 (but not its direction) we can write 

^ /(ri2,...)=n2T^/(ri2,...)- (4.3) 



dfi2 dri2 
It follows that 

-> d d d d 

di2 ■ ^23/(^2, r23, ...) = ■ T^f{ri2, r23, ■ ■ ■) = ru ■ hsT, — ^ — 7(^2, ^^23, • • •)• (4-4) 

or 12 or 23 or 12 or 23 
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Unlike the 2-body case, where the analog of Eq. ( [4.4|) simplifies because fi2 ■ r2i = —1, the 
angular factor fi2 ■ in Eq. ( ^.4| ) is not a constant, and does not average to zero when 
integrated over a sphere. Moreover, the expression in Eq. ( [4.4| ) cannot be calculated by 
separately evaluating the angular factor and the function which it multiplies, since these 
contributions are not independent. 

To further explore the angular factors we consider {fij ■ rim)R-, where the notation (■ ■ 
denotes the average over a spherical volume of radius i?, and i,j,l,m = 1, . . . , 4. {r^ ■ r34)i? 
and (fi2 ■ 'r23)ij have been evaluated numerically by randomly generating 10^ sets of 



vectors fi, r2, ra, and in a sphere of radius R = 0.5. We find 
/- - \ /(^i-^2) (r3-r4)\ 

{ri2 ■ r^lR = ( ^ ^ ■ ^ ^ ) - 0, (4.5a) 

\lr1-r2l |r3-r4|/^ 

(,-,,. ^,3), = / (4:15). fijl^)) =-0.442. (4.5b) 

The result in Eq. (|4.5a| ), which is zero to within the expected statistical fiuctuations, can be 



understood as follows: Given a sufficient number of trials (i.e., sets of vectors) we would find 
that for each configuration of 4 particles specified by fi, r2, f^, and r4 there will be another 
in which the coordinates of particles 1 and 2 are same but those of particles 3 and 4 are 
interchanged. The contributions from these two configurations to {fi2-fs4)ii evidently cancel, 
which accounts for Eq. ( |4.5a| ). By contrast, the result in Eq. ( |4.5bD is less obvious, but the 



fact that (fi2 ■ r23)_R is nonzero can be understood as follows: To start with the preceding 
argument cannot be applied to (fi2 ■ r23) r since changing the coordinates of particle 2 affects 
both fi2 and f23. That such a scalar product is non-zero when averaged over the sphere 
can be understood if we replace the unit vectors by the corresponding dimensional vectors, 
fij fij, etc. Then 

{fl2 ■ r23)/j = ((fi - f2) ■ (f2 - f3))R = (fi ■ f2 - fi ■ f3 + f2 ■ f3 - rl) R. (4.6) 



The fi ■ fj terms in Eq. (^4.6|) average to zero by an argument similar to that for (fi2 ■ 'r34)_R: 
Given a sufficient number of trials, then for each random pair of vectors f and fj, there 
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will be another in which r*j will be the same, but the coordinates of rj will be inverted 
{fj —^j). The sums of these two contributions evidently cancel, which leaves 

ri2 ■ r,,)n = {-tDr = J = --R^ (4.7) 

We note that since r2 is the coordinate of particle 2 measured from the center of the sphere, 
the upper limit on the r2 integration is R and not 2R, as in Eq. ( [4.2|) . (In the 2-body case 
considered in the previous subsection r in Eq. (|4.2| ) denotes the distance between two points 
whose maximum value is 2R.) The analytic result in Eq. ( |4.71) has been verified numerically, 
and was used ClS db check on other numerical results. 

We conclude from the above that at least some of the angular-dependent factors in V^''^ 
for k > 4 are non-zero. One can show, moreover, that even an angular factor such as 
"f^u ■ "^34; which averages to zero over the sphere when considered by itself, can give a non- 
zero contribution when multiplied by the functions which arise from differentiating Sk and 
Pk in Eq. (|3.50|) . To see this consider the term in Eq. (|3.5CI|) proportional to fi2 ■ 1^34: 

1 



ri2r23?^34r4i(ri2 + r23 + r^^ + ^41)' 



^45"! dri2 9r34 

= ri2 ■ r34 

^ 4(^12 + r|4) + + + 5(ri2r23 + r23r34 + rs^r^i + r4iri2) + 20ri2r34 + 2r23r4i 

'^12'^23'^34'^4l('^12 + ^23 + ^34 + r^l)^ 

(4.8) 

Returning to the discussion following Eq. ( [4.5b| ) we note that the contributions from two 
configurations which differ by the interchange of the coordinates of particles 3 and 4 no longer 
cancel in Eq. (|4.8| ). This is because such an interchange effectively replaces the right-hand 
side of Eq. (|;8D by 

4(r^2 + ^34) + ^24 + ^ii + 5(ri2r24 + r24r34 + r34ri3 + r^ru) + 20ri2r34 + 2r24ri3 



- ri2 ■ r34 



^12^24^^34^13 (^^12 + ^24 + ^34 + ^^13^ 



5 

(4.9) 



With the coordinates fi and r2 remaining fixed, the expressions multiplying fi2 ■ f34 in 
Eqs. (|4.8| ) and ([4.9|) are not equal in general. (The exception would be a 4-body configuration 
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in which the coordinates ri, 7^2, "^3, and formed a regular tetrahedron. However, in 3- 
dimensional space such a configuration is not possible for more than 4-bodies.) The above 
argument has been verified numerically by a Monte Carlo simulation, as we discuss in more 
detail elsewhere [ p5| . 

It follows from the preceding discussion that all of the terms in Eq. ( |3.5U| ) are expected 
to contribute to U^^\ and this has been verified numerically by means of a Monte Carlo 
simulation [^. Briefly, the symbolic program Mathematica ||36| was used at the outset 



to explicitly evaluate all the derivatives appearing in Eq. ( p. 501) , as was done in Eq. ( |4.8| 



and (|4.9| ). The resulting expression for V^^\f 12, 1^23,^34,^^1) was then evaluated for each of 
10^ configurations, where a configuration is obtained by randomly generating a set of four 
3- vectors ri, r2, ''^3, and r4. If V"/'''*(r{2) "^23; "^34; ^41) denotes the value of V"*^^^ obtained from 
the ith configuration, then 

1 Nr 

U^"^ = E ^/'V"'i2, rl3, rl4, rk) (4.10) 



where A^,- = 10^ in our calculations. We find numerically 



f/W = Af^£^y(7.7). (4.11) 

There are several implications of the 4-body results which will be useful in the ensuing 
discussion. To start with there is no theoretical reason why U^^^ should vanish, and there is 
no suggestion that it does either analytically or numerically. Secondly, the computational 
effort required to evaluate U^^^ is sufficiently great as to question the feasibility of a calcu- 
lation of even ?7(^) and U^^\ let alone U^^^ for k = C(10^^). We note that a calculation of 
U^''^ starting from a generalized version of Eq. ( |2.4| ) would require a 3fc-dimensional integral, 
and hence this approach is impractical for large k. Fortunately it turns out that for present 
purposes it suffices to determine the 0-derivative contributions Uq'^ to U^^\ for which we 
can obtain an approximate analytic result as we describe below. 

Having evaluated f/*^^^ directly by a Monte-Carlo simulation, we ask whether it is possible 
to generalize the probability density function V{r) in Eq. ( p.l2|) to the 4-body case. If this 
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were possible it would allow the Monte-Carlo evaluation to be carried out analytically, just 
as in the 2-body case. A major obstacle that must be confronted in any attempt to generalize 
V{r) is that the appropriate A;-body probability density function must depend not only on 
the separations Vij = \ri — fj\ of i and j, but also on angular factors such as fij ■ rim which 
arise from the derivative terms. Even for the 0-derivative term the generalization of V{r) to 
the k-hodj case is not known at present, as we discuss in Appendix C. For this reason, we 
develop in Appendix C the "mean value approximation" which allows the k-hodj integrals 
to be done analytically, and which we use for A; > 6. 



C. The A;-Body Potential for fc > 6 

As we have noted previously, W = J2k W^''^ = Efc U'^''^ ( fc) in Eq. (|BT]) is dominated 
by terms with k c::^ N = (9(10^^). For values of k this large, numerical evaluation of U^^^ 
is not possible, even for the 0-derivative terms. However, there is a useful bound on the 
0-derivative contribution Uq''^ to U^^^ which we derive in this subsection which, along with 
the "mean value approximation" in Appendix C, forms the basis for the ensuing discussion. 
Before doing so we present the numerical result for [/q^^ for illustrative purposes, 

Cq = 0.268(34). (4.13) 

We proceed to derive the previously-discussed bound on Uq'^ Let 'P*-^''(ri2, r23, . . . ,^"^1) 
denote the probability density for the separation of particles 1 and 2 to be in the interval 
between ri2 and ri2 -|- dri2, etc. Thus, the differential probability d'^H^''^ of a configuration 
specified by ri, r2, . . . , rfc is 

^k^ik) ^ rik)^r,2, r23, . . . , ni)drudr23 ■ ■ ■ dm, (4.14) 

where 



37 



/ dru--- drkiV^'\n2,r2s,...,rki) = l. (4.15) 
Jo Jo 

We know that the function V^''\ri2, r23, . . . , r^i) exists because it can be determined numer- 
ically by randomly generating sets of vectors ri, r2, . . . , rfc in a sphere, and then computing 
the probability d'^H^''^ of any configuration, as in Eq. (|4.14|) . In fact the first non-trivial 
generalization of V{r) in Eq. ( |2.12D , namely 'P*-^''(ri2, r23, rai), has been inferred in exactly 
this way (As noted in Sec. Ill, the fact that there is no 3-body contribution to W has 
to do with the functional form of V^*^^-*, not 'P*^^-*(ri2, r23, rai).) Combining Eqs. (|4.14| ) and 
(|3|) we find: 

U^'^ = ^ ' i'k\ f^"" dr,2 ■ ■ ■ dr,,V^'^ (r,2, r23, • • • , r,,) 



X 



(4.16) 



712^^23 ■ ■ ■ rki{ri2 + r23 H h rfci)'=+\ 

We observe that the function multiplying 'p(^)(- ■ ■) is non-negative, and is a monotonically 
decreasing function of r^j in the interval [0, 2B\. It follows that the minimum value of jf/g'^^l 
is achieved when = 2R for all Thus, 







r2R 

/ rfrfciP('=)(ri2,r23,...,r,i). (4.17) 
Jo 



Using the normalization condition, Eq. ( [4.15| ), we have 



The result in Eq. ([4.18|) represents the 0{Gp) contribution from the standard diagram 



(with both senses of the loop momentum included) as can be seen from the starting point 
in Eq. ( p.55|) . As noted in Section III and in Fig. ^, there are in 0{Gp) a total of {k — l)!/2 
topologically distinct pairs of diagrams, with each pair representing the sum of the contri- 
butions from both senses of the loop momentum. The sum of all these {k — l)!/2 diagrams 
replaces the expression in square brackets in Eq. ( ^4.16| ) by the appropriately symmetrized 
[in the k{k — l)/2 variables Vij] generalization of the standard contribution exhibited there. 
Since V^''\ri2, r23, . . . , r^i) is also a symmetric function of Vij, we conclude that each of the 
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{k — l)!/2 contributions to f/g^^ will be equal in magnitude to that arising from Eq. ( 4.16 ). 
It follows that the complete contribution to Uq'^ is simply {k — l)!/2 times that given in 



Eq. ( CT) 



It is instructive to consider the expression for ulf'^ in greater detail to understand the 
origin of the factor (fc — l)!/2. For brevity define 

u = ri2 V = r23 w = ru x = r^i y = z = ri3. (4.20) 

(4) 



From Eq. ( [4.16[ ) the complete expression for Uq is then given by 



ulf^ = — ( ^ 4! / dudv dw dx dy dzV'''^\u,v,w, x,y, z) 

^ (4.21) 



1 1 1 



uvwx{u + V + w + xY^ vxyz{v + X + y + zY uwyz{u + w + y + zY 
where the three terms in square brackets arise from diagrams (a), (b), and (c) respectively 
in Fig. 0. We note that V'^'^\u^v^w,x,y^ z) is necessarily a symmetric function of its argu- 
ments, and that the expression in square brackets is also a symmetric function of the same 
arguments. It then follows that whatever the explicit functional form of V^^'' (u, v, w, x, y, z), 
each of the 3 terms in square brackets contributes equally to the integral. For the /c-body 
contribution the expression in square brackets would contain {k — l)!/2 terms, which is the 
origin of the factor k\/2k in Eq. ( |4.19| ). 

As we have noted previously the many body contributions are dominated by terms with 
~ A^, and for these terms the bound in Eq. ( [4.19|) is useful. However, if we wish to 



obtain a closed-form expression for the sum, then the bound in Eq. (|4. 19| ) cannot be taken 
over directly, since the individual terms in the series alternate in sign, as can be seen from 
Eq. ( ^.16| ). We thus wish to replace the bound by a more precise estimate of the integral, and 
this can be done using the "mean value approximation," as discussed in Appendix C. From 
the preceding discussion, Uq'^ can be obtained by multiplying the expression in Eq. ( |U15[ ) 
by k\/2k which gives, for even k, 
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Eq. (|4.22|) is the starting point of our discussion in the next section of the combinatoric 



factors arising in a system of particles. 

V. COMBINATORICS FOR MANY-BODY SYSTEMS AND NUMERICAL 

RESULTS 

A. Combinatorics 

The expression for Uq''* in Eq. ( f4.22| ) represents the fc-body contribution to the neutrino- 
exchange energy in the approximation of retaining only the 0-derivative terms. As we noted 
in Sec. II, for a system of particles there are (^^^ such (identical) terms, and hence the 
total neutrino-exchange energy for a spherical distribution of A^ particles is given by 



even 



^ 7cR \27rV2Ry k^+^\k - 



(5.1) 



fc^4 
even 



U^'^^ is the 2-body contribution and, since it can be evaluated exactly, the full expression in 
Eq. ([4.2|) will be used. When calculating the neutrino-exchange energy in a white dwarf or 
a neutron star, 

H/(2) = f/(2) (^^^ is completely negligible and hence it can be dropped from 
the sum over k when convenient. 

The sum over k in Eq. ( |5.1| ) can be evaluated in closed form by making use of the Stirling 
approximation for k\, which is valid for large k: 



Combining Eqs. (SJ.) and ( p.2| ) we have 
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where Ine = 1. The remaining k\ will be replaced shortly by an integral representation 
which, when evaluated by the saddle-point method, amounts to a second application of the 
Stirling approximation. The sum over k can be further simplified by noting that were we to 
apply the Stirling approximation again at this stage we could write 

fcl ^fc+l/2 

^—^e''' = V2^exp[{k + l/2)\nk-k-(3/2)\nk]. (5.4) 



Since we are interested in evaluating the sum for k < N = (9(10^^), the term (3/2)lnfc 
could be dropped relative to A; In A; in Eq. (|5.4|) , and hence the factor k^^'^ in Eq. ( f>.3fj could 
also be dropped. However, a better approximation is to simply replace k^^"^ by N^^"^, noting 
that W is dominated by terms with k c:^ N. (This approximation also gives the smallest 
estimate for W.) Although the factor /c^/^ is negligible from a quantitative point of view, it 
can (with some effort) be reinstated if necessary, as we demonstrate in Appendix D. In the 
approximation of replacing k^^'^ by N^^'^, W becomes 

even 

where 



R 



(5.6) 



27r V2ei?2 ' 

Although the lower limit on the sum in Eq. (|5.5| ) is A; = 4, the summation can be extended 
down to /c = by subtracting the k = and k = 2 contributions at the end. We thus 
consider the sums, 

N 



(5.7) 



E^^^ = E^!(±^rH)^(^), (5.8) 

where in X)*'^'' the summation extends over all k. The quantity we want is J2^^\ and from 
Eqs. ( |5.7| ) and ( |5.8| ) it is given by 
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(+) , v-(-)l 



(5.9) 



To evaluate Xl''^^ we introduce into Eq. ( |5.8| ) the familiar integral representation for k\, 



k\ = / ciwe-V, 
Jo 

which gives 

The sum of the binomial series in Eq. (|5.11| ) can be expressed in closed form |37 

{a + xf, 



(5.10) 



(5.11) 



^ /iv^ 



k=0 



k 



(5.12) 



where a = 1 in the present case. Combining Eqs. ( |5.11|) and ( |5.12| ) gives 

^(±) 



due-''{l±tuTRy\ 



(5.13) 



The integral in Eq. ( [5. 131 ) can be readily evaluated by the saddle-point method We 
introduce the variable t defined by 



1 ± iuTn = ±iNt, 



(5.14) 



where the upper (lower) sign applies to ^)- can then be written in the form 

=9{N) rfte^^W, (5.15) 



where 



^(A^) = — exp 

i R 



fit) 



-— + Nln{iN) 



r 



Int 



R 



(5.16) 



(5.17) 



and tmin (tmax) correspouds to u = {u = oo) in Eq. ( |5.14|) . Using Ref. |3^, J^^'^^ is then 
given by 
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(+) v^^(Ar)e^/(*o) 



(5.18) 



- \Nf"{toW/^ ' 

where to = Tji is the saddle point, and the primes indicate differentiation. Combining 
Eqs. (I^)-(I^), we find 



^ V2^exp{N[\n{TRN) - 1] +i(iV7r/2 - 1/Tr)} 



(5.19) 



Since ^ = {Yl^'^^)*, as can be seen from Eq. (|5.11| ), J^^^'' in Eqs. ( ^.7|) and (|5.9|) is given 



by 



E 



-lf^^V2^{TRN/e)^ cos(l/rK). 



(5.20) 



Noting that only even values of contribute to X^fc; it follows that the factor (—1)^/^ is 
always real and alternates in sign. We also note that because the expression on the right 
hand side of Eq. ( |5.18| ) is the origin of both the Stirling approximation for k\ and the result 
for above, it follows that both factors of k\ in Eq. ( p.lj ) have been treated in a consistent 
way. Combining Eqs. ( |5.5| ), ( p.T] ), and ( |5.20D the neutrino-exchange energy W can be written 
in the form 



'2 2 1 



(5.21) 



where Yl^^^ and X)*^^"* are the k = and k = 2 contributions to the sum over k in Eq. 
which must be subtracted out. From Eq. (|5.7|), 



^(0) ^ ^(2) ^ ^ _ _ ^^^2 



(5.22) 



and hence 



W ~ (-1)^/2 



RN 



N 



cos(i/rij,) 



+ 



iV3 Vi? 



2 \ ^ A^(A^ 



1 



2e2 



(5.23) 
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For purposes of setting a bound on the neutrino mass the expression in square brackets is 
neghgible and can be dropped. Using Eq. ( |5.6| ) we can then write Eq. ( [5.23D as 



W ^ (-1)^/2^ { -^^P^\\os{1/Vr). (5.24) 

Although the expression in Eq. (|5.24| ) is the resuh of summing over all < in Eq. ( ^.5| ), 
it can be shown that this expression is effectively the contribution of a few terms which 
dominate the sum over k. One way of seeing this is to consider the ratio of successive terms 
in Eq. (|5.5|). Dropping the factor k~^/'^ as discussed above we find 

Tji{N - k){N - k~l). (5.25) 



As we demonstrate in Eq. ( 5.3ti| ) below, the product TjiN is of order 10^^, and hence for 



small k the ratio in Eq. ( ^.25|) is much larger than unity. It follows that successive terms in 
J2k initially make increasingly larger contributions. As noted in Sec. II, however, the rate 
at which this ratio increases is itself a decreasing function for increasing k, and eventually a 
value kmax is reached for which the ratio is unity. From Eq. ( |5.25|) kmax is given by 



r?j(Ar-fc^,,)2~l, (5.26) 

k^ax ^N- 1/Tr. (5.27) 

As noted in Sec. II, the fact that the individual terms in J2k reach a maximum is a con- 
sequence of the behavior of (^^^ for large k, and of Eq. ( |2.9| ) in particular. To calculate 
yyimax)^ the valuc of VT^^^ Corresponding to kmax, we can set kmax — N and (l/Tji) ~ since 
A^ ^ I/Tr (see below). From Eq. ( |5.5| ) we then have 

^(max) ^ (_l)iV/2^ , (5.28) 



RN \ e 

where the Stirling approximation for A^! has been used. Eq. ( p.28| ) agrees with the result 



obtained previously in Eq. ( |5.24| ), up to the factor cos(l/r/j) which is the "memory" of the 
sum over k as we discuss below. It follows from Eq. ( p. 25] ) that for practical purposes one 
can approximate W reasonably well by the single term W^"^°'^\ 
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B. Numerical Results 



We proceed to evaluate W numerically for a typical neutron star which we take to be 
the observed pulsar in the Hulse- Taylor binary system PSR 1913+16 |^D|-IT . The mass Mi 
of this pulsar is accurately known 

Ml = 1.4411(7)M0, (5.29) 

and hence the mass of a typical neutron star will be taken to be 

M = IAMq = 2.8 X 10^^ g. (5.30) 

To calculate the number of neutrons we ignore the contribution to M from gravitational 
binding energy [see Eq. (|5l3|) below], and assume that the neutron star is composed exclu- 
sively of neutrons. Using Eq. (|5.30|) then leads to 

iV = 1.7 X 10^^. (5.31) 

The radius R of the neutron star, although not directly observable, can be inferred in various 
models. We assume the nominal value R = 10 km = Rio which corresponds to a mass density 
Pm and a number density p given by 

p„ = 6.7 X 10^^ gcm-^ (5.32) 

p = 4.0 X 10^^ cm-l (5.33) 

In what follows we will assume these values of R and p, which are typical of the results 
that arise in existing models of neutrons stars [^,|^. Combining Eqs. ( p.30D -( p.32D with 
|a„| = 1/2 and reinstating h and c, we find 



7.6 X 10^^ (5.34) 



^10 



^ {G,/hc)\a \ ^ ^ ^ 1 

2nV2eRlQ 1.1x10^6' ^ ^ 
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A 



R 



{GF/hc)N\an\ 



1.2 X 10"|a„| = 5.8 X 10'^ 



10 



(5.36) 



10 



W = (-1)^/21^ (6.0 X 10-'') (5.8 X 101°)'-'''°" cos(l/r«) 
Rio ^ ^ ^ ' 



10(2xio''-57-io)(_i)^/2cog(l/r^) eV. 



(5.37) 



In Eq. ( |5.37| ) the three terms in the exponent of 10 arise, respectively, from (r^jA^/e)^, 1/A^, 
and 4/ic/i?io (in eV). For the ratio WlM(? we find 

W 



^_^)iV/2^Q{2xl0-^^«-57-76) cos(l/r^). 



(5.38) 



Leaving aside the cosine factor, to which we will return shortly, we see that the neutrino- 
exchange energy is significantly larger than the known mass-energy of the 1913+16 pulsar. 
For later purposes we note that in a neutron star W'^^'^ /Mc^ exceeds unity for k > 8, which 
is a relatively low order perturbation. Using Eq. ( p.lj) we find that for the 0-derivative 
contribution. 



2i'hc (8!)2 
nRio 810 



{GF/fic)\an\ 



10 J 



5.0 X 10^^ eV 



Mc2 



3.2 X 10^^ 



(5.39) 



(5.40) 



Another quantity of interest is the ratio of W to the total mass of the Universe Mu. We 
find: 



M[/ ~ 4 X 10^^ g 



(5.41) 



W 

Muc^ 



]^g{2x 10^9-57-99) 



(5.42) 



It follows from Eqs. ( |5.38| ) and ( |5.42| ) that the neutrino-exchange energy in a single neutron 
star, when calculated in the Standard Model, exceeds the total mass-energy of the Universe, 
and this may be termed the "neutrino-exchange energy- density catastrophe." In estimating 
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Mjj, we assume the average density pu of the Universe to be p(7 ~ 1 x 10~^^ g cm~^, which 
is the critical density corresponding to a Hubble constant of (80± 17) km s~^ Mpc~^ p^ , ^ 
and we have taken the radius of the Universe to be ^ 1 x 10^° ly (HSf. 



We return to discuss the sign of W in Eqs. ( p.24| ) and ( |5.38| ). As noted previously, the 



factor (—1)^/^ alternates in sign, and hence the sign of W depends on A^, as well as on 
Vr through cos(l/r/j). Either sign of W leads to a catastrophic outcome: W > Q would 
correspond to a large repulsive force against which the neutron star would be unstable, 
whereas W < Q would cause the neutron star to collapse to a black hole. In either case 
neutron stars as we know them would not exist. In this context it is helpful to recall that 
the classical gravitational potential energy Vgbav is approximately given by 

Vgrav ^ ^ , (5.43) 

where Gn = 6.67259(85) x 10^^^ m'^kg^^s^^ is the Newtonian gravitational constant. For 
the assumed values of M and in Eqs. (|5.3CI|) and ( |5.31|) we find 

I Vgrav I 



Mc2 



0.124, (5.44) 
117 MeV/nucleon. (5.45) 



N 

It follows from Eqs. (|5.38|) , (|5.44| ), and ( ^.45| ) that the neutrino-exchange energy, whatever 



its sign, would dominate over Vgrav i which is the largest contribution to the binding energy 
of a neutron star. We note in passing that since Vgrav is negative, the effect of including 
gravitational binding energy in Eq. ( ^.30| ) would be to increase N in Eq. (|5.31| ), which would 
make the neutrino-exchange energy density even larger than the value quoted in Eqs. (|5.37|) 
and (ICTp . 



We turn next to the factor cos(l/r^) = cos{2T{\/2eB? /Gpan). This is similar to the os- 
cillatory factors that arise in other many-fermion systems, an example being the Ruderman- 



Kittel interaction |4^. To better understand the significance of this factor it is instructive 
to ask what the expression for W would be if all terms in X]fc had the same sign (which we 
assume to be positive for illustrative purposes). The analog of Eq. (|5.7| ) would then be 
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(e) N /]\f\ 

S =E*n(,). (5.46) 



even 



and hence from Eqs. ( |5.7| ) and ( |5.4ti| ), 

E =E (r^) = E (r/?--^r^)- (5.47) 



It follows from Eq. ( p.47[ ) that W in Eq. ( |5.24| ) would be replaced by W where 



and hence the difference between W and is the replacement 

(-1)^/2 eos(l/r^) ^ cosh(l/rR). (5.49) 

Since cosh(l/r/j) never vanishes, W is always nonzero for the values of the parameters we 
are assuming. By contrast cos(l/rij) can be zero for some set of parameters. Hence W can 
in principle vanish, but only when 2/(7rr/j) coincides with an odd integer to ~ 10^^ decimal 
places. Not only would this be unphysical for even a single neutron star, it could hardly 
be supposed that cos(l/r/j) would vanish in this way for each of the more than 500 known 
pulsars From the recent catalog of 558 pulsars by Taylor et al. we note that there 



is a significant variation in the period P and spin-down rate P both of which depend on the 
internal structure of the pulsars. Absent an overarching (and presently unknown) symmetry 
principle which would force W to vanish, it is difficult to imagine that this could happen 
coincident ally in every case. This argument is further strengthened by noting that the 
catalog of Taylor et al. specifically excludes accretion-powered systems which are detected 
by their X-ray emissions. For these, the in-falling matter would continuously change the 
parameters of the pulsar at a level that would preclude the possibility that cos(l/r/j) would 
always be zero. 

In the preceding discussion we have considered the effects of the neutrino-exchange energy 
in an idealized non-rotating neutron star. Including the effects of rotation would further 
bolster the preceding arguments, since the variation in P implies a variation in the internal 
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structure of the pulsars, as noted above. Proceeding further one could also ask whether a 
neutron-star could even come into existence in the presence of an energy density as large 
as would arise from neutrino-exchange. Although the answer to this question could lead 
to a more stringent limit on neutrino masses than is obtained from our "static" picture, 
such dynamical considerations are beyond the scope of the present paper. However, even 
in the absence of a detailed calculation one can argue that normal stars could not evolve 
into neutron stars in the presence of neutrino-exchange. This follows by noting that as soon 
as a star evolves to the stage where (TuN/e) in Eq. ( ^.231) exceeds unity, its mass-energy 



not only becomes unphysically large, but also alternates in sign as each pair of neutrons is 
accreted. Hence if at some stage (—1)^/^ were —1, the addition of 2 neutrons would require 
overcoming an unphysically large repulsive energy barrier, which thus acts to prevent further 
accretion. 

Having demonstrated that the oscillatory factor cos(l/rij) cannot prevent the neutrino- 
exchange energy- density catastrophe, we study some of the implications which follow from 
the presence of this factor in the expression for W. To start with, the sign of the 0-derivative 
contribution to W in Eqs. ( |5.24| ) and ( |5.37| ) depends on the product (— 1)^/^ cos(l/r/j). To 



determine the sign and magnitude of cos(l/rj:j), Tji itself would have to be known to an 
unphysically large number of decimal places. (For illustrative purposes, if in Eq. (^.35|) 



were exact, then cos{l /Tji) = +0.613 [^.) Hence as a practical matter, the overall sign 
of W for a given neutron star is difficult to determine, but is also relatively unimportant 
given that either sign leads to a difficulty. There is, however, an interesting implication 
of the variation in sign of as a function of N and Tfi: Rather than limit our attention 
to the parameters of the ~ 500-600 known neutron stars that are pulsars, we can view 
any spherical volume of radius r inside a neutron star as a collection of N{r) neutrons to 
which the preceding arguments apply. In order to have a stable static neutron star all such 
subvolumes must be in equilibrium, notwithstanding the variation of the sign and magnitude 
of W{r) with r. By pursuing this argument, we are led in Section VIII to the bound on m 
quoted in Eq. (|8.12|) below. 
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We conclude the present discussion by presenting a physical interpretation of the factor 
;-l)^/2cos(l/rH). Combining Eqs. QOSl ) and j ^AEj) we can write 



iy('""^)cos(l/rjj), (5.50) 
W ~ (_i)^/2|y("^<^^) cosh(l/r^). (5.51) 
We see from Eqs. ( ^.50|) and ( ^.51| ) that the overall scale of both W and W is primarily 



determined by the sing le term VT^™"^) in Eq. (|08|) . The factor cos{1/Tr) in Eq. ( ^M) 



which acts to suppress W, can then be understood as the "memory" of the cancellations 
arising from the remaining terms in X)''^^- By way of contrast the factor cosh(l/r/{) in 
Eq. ( ^.51|) enhances the contribution from W^'^"'^\ and thus reflects the effect of coherently 



~(e) 

adding all the same-sign contributions in J2 ■ For practical purposes the magnitude of 
the neutrino-exchange energy- density is thus determined by W^'^"'^\ for which the bound in 
Eq. ( [4.19|) can be used. 



While on the subject of cancellations it is worthwhile to consider another possible source 
of cancellations, namely the effects of the Pauli exclusion principle for neutrinos. The Pauli 
effect for neutrinos can enter in two ways: Firstly, the exchanged pairs can interfere with 
one another and, secondly, the exchanged neutrinos could interfere with a possible neutrino 
sea produced by the neutron star or white dwarf. With respect to the first possibility, the 
exchanged neutrinos do in fact interfere with one another, but these effects are rigorously 
taken into account from the outset in the Schwinger-Hartle formalism, and at all stages in 
the present derivation. Since the Schwinger formula in Eq. ( P34| ) is simply a convenient 
way of summing all the one-loop Feynman amplitudes (which clearly incorporate the Pauli 
exclusion principle), all such effects are automatically included. For example, the factor of 
(—1) from the closed neutrino loop, which is a consequence of the anticommutativity of the 
neutrino field operators, is already built in. Moreover, the expression for W in Eq. (|5.24| ) 



is the result of coherently summing the fc-body contributions, with all the relevant phases 
determined by the Pauli principle. Hence all coherent effects of the Pauli exclusion principle 
for the exchanged neutrinos are accounted for in the present results. 
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Consider next the possibility that neutrino-exchange is somehow suppressed by the pres- 
ence of a sea of physical neutrinos in a neutron star. It is easy to demonstrate that such a 
suppression cannot take place, because the density of physical neutrinos in a neutron star is 



Here pnuc is the nuclear density, which is taken to be pnuc = 2.8 x 10^^ g cm~^, and pm is the 
neutron star density. For E^, = 10 eV, the energy scale set by Gppn, we find ^ 1 x 10^° 
km, which means that a neutron star is transparent to low energy neutrinos [|^. Since 
low-energy neutrinos cannot be trapped inside a neutron star, there is no possibility of a 
neutrino sea in a neutron star. The few neutrinos that are present (in transit) at any instant 
are obviously incapable of reducing the effective Fermi constant by a factor of order 10^^, 
which is what would be needed to resolve the energy-density problem. Finally we recall 
from Table ^ that neutrino-exchange leads to a large energy-density in white dwarfs as well 
as in neutron stars. Hence this difficulty cannot be resolved by invoking a physical neutrino 
sea, since the mean free path of a neutrino in a white dwarf would be even larger than in a 
neutron star, so even fewer neutrinos could be trapped to form a neutrino sea. 

Returning to Eqs. (|5.24|) and ( ^.28| ) we note that these results arise from the 0-derivative 
contribution Uq'^ in Eq. ( ^.1|) . We consider next the contributions to W from the derivative 
terms in Eqs. ( p.55D -( ^.58| ). Since these terms introduce no new dimensional factors, their 
contribution to W must have the same dependence on {Gpanl R^Y in Eq. (|0|), although 
the fc-dependent coefficients will in general be different. We consider three possibilities: a) 
The net contribution from the derivative terms is smaller than that from the 0-derivative 
terms. In this case Eq. ( [5 . 1| ) is an adequate approximation to W , and the previous conclu- 
sions apply as before, b) The net contribution from the derivative terms is larger than that 
from the 0-derivative terms, a possibility suggested by the 4-body results. In this case the 
neutrino-exchange energy- density catastrophe poses an even greater difficulty than before. 



essentially zero From Ref. [^] we note that the mean free path A„ for the scattering 



of Uf, from neutrons is given by 




(5.52) 
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and resolving the problems raised by the 0-derivative contribution is a necessary first step. 
It is evident that whatever fc-dependent coefficient replaces (A;!)^/(A;'^+^) in Eq. ( |5.1| ), the 
fc-body contribution arising from the derivative terms will necessarily be proportional to the 
same factor 



that appears in the 0-derivative contribution. Hence whatever mechanism resolves the 
neutrino-exchange energy-density catastrophe for the 0-derivative terms, will also work for 
the derivative terms, as we discuss in greater detail in Sec. VIII below, c) The last possibility 
is that the derivative terms cancel the 0-derivative contribution so as to reduce \W/Mc^\ in 
Eq. (|5.38|) to a number of order unity (or more realistically to ~ 0.1). Such a cancellation 
would have to occur in each of the 500-600 known pulsars, and (as we discuss below) in each 
of ~2000 known white dwarfs. Since the derivative terms have a different dependence on fij 
than the 0-derivative terms, they would also depend differently on the matter distribution 
in the neutron star, were we to allow for a varying density in a more realistic calculation. 
Thus, to avoid the neutrino-exchange energy- density problem would require at a minimum 
an almost exact cancellation between the derivative and 0-derivative terms in each pulsar, 
notwithstanding the fact that the former depend differently from the latter on the matter 
distribution, which itself varies from one pulsar to another. As noted previously, such can- 
cellations would be unphysical absent some (presently unknown) symmetry principle. We 
are thus led to the conclusion that the derivative terms can exacerbate the energy-density 
problem arising from the 0-derivative terms, but they cannot resolve it in the framework of 
the present (Standard Model) calculation. 

As noted in the preceding paragraph, neutrino-exchange leads to an energy-density catas- 
trophe in white dwarfs as well as in neutron stars. Although neutron stars provide the most 
stringent lower bound on m, the weaker bound from white dwarfs is nevertheless important 
because its existence makes it less likely that this problem can be resolved without intro- 
ducing massive neutrinos. For example, the fact that white dwarfs contain roughly com- 




(5.53) 
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parable numbers of electrons, protons, and neutrons, rules out any attempt to address the 
problem in neutron stars by a mechanism which somehow suppresses the neutrino-neutron 



coupling. A 1987 catalog by McCook and Scion lists 1279 white dwarfs which have been 
identified spectroscopically, and the current catalog contains more than 2000 entries. To 
demonstrate that the energy- density arising from neutrino-exchange in white dwarfs is also 
unphysically large, we consider two representative white dwarfs whose masses and radii are 
reasonably well-known |5l[]. The relevant parameters for these two white dwarfs, Sirius B 



and 40 Eri B, are summarized in Table |l[ Following Ref. we assume that the interiors of 



these white dwarfs are predominantly carbon and oxygen, in which case these stars contain 
approximately equal numbers of neutrons, protons, and electrons as noted previously. The 
combinatorics of the many-body diagrams involving three types of particles are somewhat 
complicated, but for present purposes it is sufficient to evaluate the purely electronic con- 
tribution, noting from Appendix A that electrons have the strongest coupling to neutrinos. 
The previous formalism for neutron stars can then be taken over immediately and leads to 
the results given in the last two lines of Table |I| where. 

Here Nf. is the total number of electrons, a^. = (2sin^ 6^ + 1/2) ~ 0.964 is the e-z/ coupling 
constant [see Eq. ( |A9c|) ], and we assumed for illustrative purposes that the electron density 



is spatially constant. It follows from Table |l|that for both Sirius B and 40 Eri B the energy 
arising from neutrino-exchange exceeds the mass-energy of the Universe, just as in the case 
of neutron stars. 

The results in Table raise the possibility that the existence of a large neutrino-exchange 
energy-density may be not be limited to neutron stars and white dwarfs. From Eq. ( |5.24| ) 
we see that for a single particle species j = n, e, or p, this happens when 

where aj is the coupling to neutrinos. For the Sun the contribution to from electrons 
alone slightly exceeds unity, which suggests that the very existence of the Sun, and by 
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extension life on Earth, is evidence for the necessity of massive neutrinos. Whether or not 
the Sun itself proves to be unstable against neutrino-exchange can only be determined by 
more detailed calculations. However, the possibility that other stars exist for which this is 
true indicates that the problems arising from neutrino-exchange can only be resolved by a 
universal mechanism that would apply to large numbers of rather different celestial objects. 



VI. ALTERNATIVES TO MASSIVE NEUTRINOS 

We consider in this section several alternatives to the introduction of massive neutrinos 
as a means for resolving the neutrino-exchange energy-density problem. It will be argued 
that no presently known mechanism except for massive neutrinos is compatible with existing 
experimental and theoretical constraints. 



A. Deviations of the Neutrino Couplings from the Standard Model Predictions 



As noted in the Introduction, the magnitude of W in Eq. ( p. 381 ) is determined by dimen- 



sional arguments and combinatoric considerations, and hence is insensitive to the detailed 
form of the neutrino coupling. To date there are no known discrepancies between the predic- 
tions of the Standard Model and experiment (see Appendix A), and the agreement is at the 
level of a few percent or better . It would thus be difficult to understand how the effective 
Fermi constant Gi^lanl could differ from the Standard Model value by a factor of order 10^^, 
which is what would be required to avert a large energy-density. Similar considerations 
apply to the possibility that the neutral current interaction contains small admixtures of S, 
P, and T in addition to the usual V, A currents. 

In this connection we recall that the mass density pm in Eq. ( |5.33| ), pm = 6.7 x 10 



14 



g cm , is comparable to nuclear density, pnuc = 2.8 x 10 g cm . It follows that the 



success of conventional theory in explaining weak interactions in nuclei |Q strongly suggests 
that the weak neutrino-hadron and neutrino-electron couplings in neutron stars do not differ 
significantly from the predicted Standard Model values. This is supported by the general 
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agreement between theory and observation for such processes as supernova formation ^3 . 
We further observe that even smaller renormalization effects are expected in white dwarfs, 
where a large neutrino-exchange energy-density nonetheless exists. Finally, we note that 
various astrophysical arguments support both the spatial and temporal constancy of 
53]. Taken together these arguments make it unlikely that the present difficulties can be 



resolved by modifying the neutrino coupling constants. 



B. Cancellations Among Vg, Vn, and 



The expression for W in Eq. ( ^.24| ) represents the contribution from a single neutrino 



species, z/gi/e for example. It might be argued perhaps the contributions from u^Ue, 
and VrJ^r could cancel amongst themselves in such a way as to reduce to a physically 
reasonable value. However, such a scenario is in conflict with e-/i-r universality, which 
implies that the couplings of z/g, v^-, and z/^ to n, p, or e have the same signs as well as 
the same magnitudes, and hence cannot produce the necessary cancellation. Experimental 
support for e-/i-r universality comes from a number of sources including the equality of the 
partial decay rates V{Z^ where £ = e, /i, r 

C. Breakdown of Perturbation Theory 

It might be suggested that the fact that W is unphysically large indicates that our use 
of the Schwinger formula in Eq. (P34|) is not valid for some reason. This is a possibility 
that cannot be completely excluded at present, particularly in light of recent interest in the 
breakdown of perturbation theory in high orders . We can argue, however, that this is not 
likely to lead to a resolution of the energy- density difficulties for reasons we now discuss. To 
start with, the expression for W in Eq. ( [5.24| ) is the sum of a finite number of contributions, 
each of which is well-behaved. Moreover, for > 6 the one-loop A;-body diagrams are 
finite, even in the absence of any regularization scheme. It follows that the magnitude of W 
cannot be explained as an artifact of the manner in which the potentials are extracted from 
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one-loop amplitudes. Finally, the actual k-hodj neutrino-exchange potentials which arise 
from the Schwinger-Hartle formalism are themselves well-behaved; it is only when particular 
values of the parameters Gp, N, R, p, and m are used (e.g., m = 0) that unphysical results 
emerge. This suggests that the resolution of the present difficulties lies in establishing what 
the correct value of m is rather than in finding a possible breakdown of perturbation theory. 
The previously discussed calculation of the Coulomb energy of a nucleus provides a useful 
analogy. Suppose that it was believed, for whatever reason, that for a nucleus with A 
nucleons the nuclear radius R in Eq. (|2.6|) was given hy R = A^/^o, with ro = 0.01 fm. The 
resulting Coulomb repulsion would be large enough to destabilize all nuclei, so matter as we 
know it could not exist. Although one could attempt to resolve this "Coulomb catastrophe" 
by examining higher-order effects in perturbation theory, a more reasonable approach would 
be to first ask whether another value of the parameter rg, say tq ~ 1 fm, was compatible with 
experiment. Analogously, we demonstrate in Sec. VIII below that if m ^ 0.4 eV/c^, then 
neutrino-exchange no longer leads to an energy-density catastrophe. This path has clear 
observational implications, and should these be shown to be incompatible with experiment, 
then a re-examination of perturbation theory would be appropriate along the lines we now 
consider. 

As we have already noted, the final expression for W in Eq. ( ^.24| ) is the result of 
summing the contributions from a finite number of many-body potentials, each of which 
is well-behaved. It follows that if there is indeed a breakdown of perturbation theory, it 
cannot be due simply to a failure of the perturbation series to converge. One possibility 
which has been discussed recently |^5[| is the rapid growth of the cross sections for producing 



large numbers of scalar particles. It has been conjectured that the combinatoric factor n\ 



which appears in the amplitude for producing n particles in the final state could overcome the 
smallness of the coupling constant in weakly coupled theories, and thus produce observably 
large effects at high energies. This has led to an examination of the need for unitarizing the 
tree-level amplitudes from which these effects arise, and more generally, to a discussion of 
whether such effects can in fact be seen. 
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Although the origin of the combinatoric factor n! which is responsible for the growth 
of the tree-level contributions is similar to that for the factor {k — l)!/2 in Eq. ( f4.19| ), the 



neutrino-exchange problems depend critically on the presence of the factor A^! in (^^^ , which 
has no analog in the work described in Ref. |^5[. Moreover, we have seen in Eq. ( [5.40[ ) that 



the large neutrino-exchange energy-density arises in 0{G^p), which is a relatively low order of 
perturbation theory compared to the effects considered in Ref. [^]. Finally, the self-energy 
of a neutron star arising from neutrino-exchange is a low-energy phenomenon, in contrast 



to the production processes in Ref. ||55|. Hence the unitarization mechanism considered in 
that context would not be relevant here. For all these reasons it appears unlikely that a 
breakdown of perturbation theory along the lines contemplated in Ref. could resolve the 
neutrino-exchange energy- density catastrophe in a neutron star or white dwarf. Nonetheless 
this remains an interesting avenue for future exploration. Since at present the only known 
viable mechanism is the assumption that neutrinos are massive, we explore the implications 
of m 7^ in the following sections. 

VII. MANY-BODY INTERACTIONS FROM THE EXCHANGE OF MASSIVE 

NEUTRINOS 

In this section we apply the formalism for massive neutrinos developed in Appendix E 
to the problem of calculating the many-body neutrino-exchange potentials. Our objective 
is to demonstrate that when m 7^ 0, the neutrino-exchange energy can be approximated by 
replacing 

l/R" ^ e-'"^//^^ (7.1) 



in Eq. ( [5.24| ). As noted in Appendix E this leads to "saturation" of the neutrino-exchange 



interaction, and ultimately to a bound on m as we describe below. 

When m the couplings of neutrinos to other matter may be different from the 
predictions of the Standard Model. Although the details of whatever Hamiltonian would 
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emerge cannot be fully discerned at the present time, it is sufficient for present purposes to 
consider a simple phenomeno logical model in which z/e, u^, and z/^ remain as the eigenstates 
of the Hamiltonian, with each having a nonzero mass. The modification of the neutrino 
couplings to matter can then be parameterized as in Eq. ( |7.2| ) below. 



A. The 2-Body Potential when m ^ 

The expression for W^'^^ when m ^ can be obtained from the m = expression in 
Eq. (|3.3| ) by replacing S p\f 12, E) with the expression given in Eq. ( |E18|) . In addition we 



allow for the possibility that the neutrino coupling may not be pure V-A when m 7^ by 
substituting 

(1 + 75)^(1 + 675), (7.2) 

where 6 is a constant (which may be different for u^, ly^, and u^-). If there is also a modification 
of the overall strength of neutrino coupling, this can be accommodated by appropriately 
redefining Gp- After calculating the trace of the Dirac 7-matrices, the analog of Eq. (p.lO|) 
can be written in the form 

^(2) = I (^^^ ' I d'x, I p2 d'x2 dE{il + b')E' - [(1 + b')d,2 ■ 821 - (1 - 6^)m2]} 
2 



I V 1 



X — exp[i(ri2 + r2i)V^2_^2]_ (73) 

KA-Kj ri2r2i 

Following the discussion in Sec. Ill, and using Eq. ( [E20| ), the 2-body potential V"^^) IS given 
by 

lOvT'^ J ri2r2i 

where z = ri2 + r2i- The expression in curly brackets in Eq. ( [7.4[ ) can be simplified by using 
Eq. ( ^.171 ) and the recurrence relation in Eq. (^3Up which leads to 



lOvr^ I ri2r2i 9ri2 9r2i 
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ri2r2i 



2m' k . (7.5) 

ri2?^2i 



After the derivatives in Eq. (|7. 5|) are explicitly evaluated, and use is made of Eq. ( [E33| ) , we 
find 



^irf)(2mr)-^iri(2mr) 



nr 

—Ki[2mr) 



(7.6) 



where r = ri2 = r2i. Note that the contributions proportional to f^j , which arise from 
the derivatives and from Fq^\ cancel against each other. Using Eq. ( [E35D Ki^\2mr) can be 
eliminated in favor of Ko{2mr) and Ki{2mr), and this leads to the final expression for 
when m 7^ 0: 



47r3 



TTl TTl 

—Ki{2mr) + —Ko{2mr) 



+ —Ki{2mr) 



(7.7) 



The result in Eq. ( [r.7|) , which is exact, can be compared to the corresponding result in 
Eq. ( |3.19|) for the massless case by expanding Ki{x) and -^'0(2;) for x ~ 0: 



Kq{x) ~ — Inx, 
Ki{x) ~ 1/x. 



(7.8) 
(7.9) 



The only contribution which survives as m — > is from the term proportional to and 
this gives 



GUI (1 + 



(7.10) 



When 6=1, which is the value appropriate to the massless case, Eq. ( [7.10|) reproduces 
Eq. (p. 19 ) as expected. As noted in Appendix E, however, the regime of interest here is 
when 2mr ^ 1, in which case the asymptotic expression in Eq. ([E27|) for K^{2mr) leads to 
the anticipated exponential fall-off of V^'^^r). 



B. Absence of Odd-A; Contributions when m ^ 

In this subsection we demonstrate that for m 7^ it remains the case that there are no 
contributions to W for odd k. Using Eq. ( |3.29aD we have 
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C "^(7 ■ 7] — m)C = (—7 ■ 7] — raf 



(7.11) 



from which it follows that 



C ^'S'^m(^u) E)C — Spl^ {—nj, —E) — Spin {^jii —E). 



(0)T, 



(7.12) 



This is the same relation that holds in the massless case [see Eq. ( p.3(]| )], and because the 
sign of E is immaterial for the reasons discussed previously, we can write as before, 



C ^ Spl^{riij)C — Spli (?7ji)- 



(0)T, 



Under C the modified neutrino coupling transforms as 



(7.13) 



c-\i + irf,)c = {i + in,y 



(7.14) 



and since both the coupling and the neutrino propagator behave under C just as they do 
in the massless case, the previous conclusion that there are no odd-k contributions follows 
when m 7^ as well. 

It is instructive to illustrate the preceding conclusion by considering the 5-body contri- 
bution as an example. The contribution from the standard diagram is given in the static 
limit by 



2^ (^) /l^^Tr{7,Xl+&75)41(51)7.(l+&75)41(45)7p(l + &75)41(34) 
X 7a(1 + 675)S'^°i(23)7,,(l + b'y5)Spl,{12)^Tf,^pXu{xi, X2, X3, X4, X5). (7.15) 



In the static limit the expression in curly brackets is given by 



tr [Eq. (7.15)] = tr [(1 - 675)4°i(51)(l + &75)5Z(45)(1 - 675)^^(34) 
X (1 + 675)41(23)(1 - 675)744°i(12)] , 

where Spl^{51) = 74^^(51)74. Combining Eqs (|7.13|) -( [7.16| ), the trace in Eq. ( [7.16|) can 
be re-expressed in the form 



(0) 



:(o) 



(7.16) 



tr[Eq. (7.15)] = -tr (1 - 675)41(-12)74(1 - &75)4I(-23)(1 + 675)4I(-34) 



(0) 



X (1 - 675)41(-45)(1 + 675)4I(-51) 



(0) 



(7.17) 
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The contribution from the diagram with the reversed loop momentum can be written down 
in a similar manner, and the analog of the right-hand side Eq. ( [7.17| ) is 

tr [(1 + 675)41(21)74(1 + 675)41(32)(1 - b^,)sP^{A3){l + b^,)SP^{5A){l - 675)41(15); 

(7.18) 

:(o) f - c(o) 



Since Spl^{—ij) = Spl^{ji) it follows that when the contributions in Eqs. ( |7.17| ) and ( |7.18| ) 



are added only those terms containing odd powers of 6 survive. All such terms can be 
eventually reduced to the trace of a product of 7-matrices containing a single 75, and from 
Eq. (|A3|) it is seen that any such trace is always proportional to e^^^xp- By the arguments 
given previously in the massless case, all terms proportional to e^,yXp average to zero when 
integrated over a spherically symmetric matter distribution. This leads to the conclusion 
that there is no 5-body contribution to W, and by extension no contribution for any odd k, 
even when m 7^ 0. 

C. The 4-Body Potential for m / 

In this subsection we obtain the explicit form of the 4-body potential, which is then used 
to infer the dependence of the general A;-body potential on the neutrino mass m. This result 
forms the basis for the bound on m that we derive in the following section. 

As in the 2-body case we begin with Eq. ( |B34| ) and substitute (1 + 675) for (1 + 75), and 
the massive neutrino propagator in Eq. ( [E18D for the massless one. The formalism of Sec. Ill 
can then be taken over directly and, after the Dirac traces are evaluated, the expression for 
the 0-derivative contribution from Fig. ^(a) to the 4-body potential is given by 

Vo^'\fu, f23, f34, f4i) = 8 dE [E'cib) + m'E'c2ib) + m%{b) 

X A^^(fi2, E)Ai.„(f23, E)Ap^(r34, E)/\Fm{Ul. E), (7.19) 



where 



Apmiru, E) = exp (ir^2^ E^ - mA , (7.20) 

47rri2 ^ ^ 
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c,(b) = {l + b'y + Ab' ; C2{b) = 2{l-b'){3 + b' 



co{b) 



'1-h 



2\2 



(7.21) 



The expression for V^^^^ given in Eqs. ( [7.19| ) — ( [7.21[ ) is the sum of the contributions from 
the standard diagram with both senses of the neutrino loop momentum. This introduces a 
factor of 2 in Eq. ( [7. 191 ), and combined with a factor of 4 from the Dirac trace, accounts 
for the factor of 8. As will be clear from the ensuing discussion, it is sufficient for present 
purposes to evaluate the 0-derivative contribution to V^^\ from which the dependence of 
on m can be inferred. From Eqs. ( [7.19|) and (|7.20|) , 



Apmiru, E) A Fm{r 23, E)AFm{rM, E)AFm{ni, E) 



-V 



All J ri2r23r34r4i 
X exp 



_ m^{ri2 + ^23 + ^34 + r4i) 



(7.22) 



which allows Vq^'^ to be expressed in the form 

1 



(4) 



TT V V2 



47r/ ri2r23r34r4i 



Ci{h)F^{z) + m'c2{b)F2{z) + m%{b)Fo{z) 



(7.23) 



Here z = (ri2 + r'23 + r34 + r4i), and the functions Fn{z) are defined in Eq. ([E20| ). Combining 
Eq. {^2^ with Eqs. ( [E30|) and (|E3g) then leads to 

w _ _,4 /G^y /^y — I — 

TT V V2 / V47r/ ri2r23r34r4i 
X {c,ib)Fl^'\z) - [2c4(6) + C2ib)]m'Fl,'\z) + [04(6) + 02(6) + Co(6)]m^Fo(^)} . (7.24) 

It is worth noting that there are mass-dependent contributions even in the limit of a pure 
V-A coupling (6 = 1). This follows from the observation that although 02(1) = Co(l) = 0, 
each of the mass terms receives a contribution from 04(1) = 8. Combining Eqs. ( [7.24|) and 
( PDSOp the complete expression for can be written as 



^0^ 



(4) 



.4 / Gi^a^X^ / ? ^ ^ 



7T \ ^/2 

X < €4(6) 2i 



-Y 



1 



4nJ ri2r23r34r4i 
2m* 



+ 



12m2\ ^ ^ / 5 7m3 4!\ ^ , , 
— — Ko{mz) + m'^ H — + m— Ki{mz) 
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- [2c4(6) + C2(6)] 2z 



. / , ,2! 



Ko{mz) + [m!^ + m^—\ Ki{mz) 



z \ z'^ 



+ [C4(&) + C2(6) + Co(6)] 2inV>Ki{mz)] . (7.25) 

As in the 2-body case we wish to check V^*-^^ in the hmit m — > 0. Using Eqs. ( [7.8| ) and 
( [7. 91 ) we note that since x'^lnx ^ as x — all the terms containing KQ{mz) vanish in 
the m = limit. Among the terms containing Ki{mz) only the term proportional to l/z"^ 
survives, and when 6=1 this gives 

Vo''' ^ - ( ^ ^, (7.26) 

which agrees with the 0-derivative contribution in the massless case as given in Eq. (p.50|) . 



When m 7^ the spatial dependence of Vq^^ receives contributions from terms propor- 
tional to m^, m'^/z, m? / z^ , and m/ z^, each multiplying either KQ{mz) or Ki{mz). We note 
from Eq. ( [E27| ) that when mz ^ 1, K,y{mz) can be approximated by 



(4) 



SO that in the asymptotic regime Vq contains terms of the form 

nmz ^, / m? m? m 1 



e-™M— . (7.28) 
2 \ z z^ z-^ z^ z^ J 

For values of mz where the exponential makes a significant contribution to W one can 



approximate \J-Kmz/2 by unity. It follows that for the term proportional to l/z^, which is 
the origin of the m = result in Eq. ( [7.6|) , the most important consequence of a nonzero 
neutrino mass is that 

-1 -mz 

- ^ (7-29) 

z^ z° 

as noted in Appendix E. Since z = (ri2 + r23 + r34 + r^i) the effect of the mean value 
approximation is to replace z by AR so that 



1 /e™^\ 1 1 / e~"^^\^ 



ri2r23r34r4i \ z^ J R \ 
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(7.30) 



We conclude from Eq. (|7.3CI| ) that for the contribution to W^'^^ arising from the primary 
effect of a nonzero neutrino mass is that l/R"^ is replaced by exp(— mi?)/i?^ as expected. 

It is straightforward to show that this result can be generalized to the A;-body case. 
In order Gp the contribution which reproduces the m = result arises from the terms 
proportional to E'^, which leads to the function Fk{z) in Eq. ( [E20| ) . The generalizations of 
the recurrence relations in Eqs. ( [E30| ) and ( [E32D eventually express Fk{z) in terms of the 
k-th derivative K'^\mz) of Ki{mz) by using Eq. ( |E33| ). These derivatives can be evaluated 
from Eq. ( [E35| ), 

K^l\z) = -Ko{z) --^K,{z) (7.31) 

by repeatedly using Eq. ( [E34| ) . Among the terms that contribute to K^^ will be one which 
arises from successive differentiations of I/2, and this produces a term proportional to 
{—l)^k\/z^. When all the appropriate factors are included we find 

E'' ^ Fk{z) {-ifFil'\z) = {-i)Him''+^ K^^\mz) 

= — —-2imzKi(mz) + ■ ■ ■ 



where the dots indicate the remaining contributions to Ki'\mz). The coefficient of the 
expression in parentheses is the massless result in Eq. ( p.l4| ), while the remaining factor is 



the leading m 7^ modification of the massless result in the asymptotic regime. It follows 
that in the mean value approximation, where z ~ kR, the modification of the massless result 
can be approximated by the factor 



-^e"™^ ^ W — - — e""''-^. (7.33) 

When calculating W the contribution from the exponential is dominated by values of the 
argument near unity, in which case (7rm/ci?/2)^''^ is also of order unity and can be dropped. 

(k) 

We conclude that when m 7^ the contribution to Uq from the term we are considering in 
Eq. ( [4.22|) is approximately given by 
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The contributions from the other z-dependent terms in Eq. ( |7.25| ) can be treated in 
a similar manner. Since each term contains either KqIitlz) or Ki{mz), both of which are 
proportional to e""*^, it follows that each contribution to Uq^'' will contain the damping factor 
exp{—mkR) which leads to saturation of the neutrino-exchange forces. The contributions 

(k) 

to Uq from this factor are of the form 



where a and (3 are integers, and where Ck is a /c-dependent coefficient which can in principle 
be determined for each such term. In practice this would be not only tedious but also 
unnecessary, since any combination of such terms leads to roughly the same limit on m. For 
purposes of deriving the limit on m it is helpful to note that for each k there must be at 
least one term corresponding to a = (3 = 0, since this is the only term which reproduces the 

(k) 

known contribution to Uq in the m = limit. In the next section we use the preceding 
results to derive a lower bound on the mass of neutrinos. 



VIII. BOUND ON THE NEUTRINO MASS 

In this section we derive the actual bound on m which follows from the assumption that 

the mechanism for resolving the neutrino-exchange energy-density catastrophe is a nonzero 

value of m. A simple bound on the mass m of any neutrino can be inferred from the 

observation that if the Compton wavelength of the neutrino were larger than the radius 

RiQ of the neutron star, then neutrino-exchange forces would behave as if neutrinos were 

massless. This leads to the estimate, 

Ti Tic 

— <Rio mc^ > — = 2x10-^^ eV. (8.1) 
mc Rio 

However, a mass this small would be insufficient to prevent a smaller subvolume of the 
neutron star from having an unphysically large energy density. Specifically, if the neutrino- 
exchange energy in a given subvolume were to exceed the known mass (and by extension all 
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other interaction energies) contained in that subvolume, then the subvolume would become 
unstable against small perturbations, in analogy to Earnshaw's theorem for electrostatics 
56| , |57[| . Following Ref. |Q, let $(xo) denote the potential energy of a neutron located in 



the neutron star at a point xq, at which no other sources are present. If the neutron is 
displaced infinitesimally from xq to Xq + 6x, where 6x = {Sxi, 6x2, Sx^), then the change in 
its potential energy is given by 

A<l> = $(fo + Sx) - ^(xo) ~ 6x ■ V$(fo) + i E Sxi6xj^^^^. (8.2) 

i,j=i ' J 

For the neutron to be in equilibrium to start with, it must be the case that the net field 
at Xq must vanish, and since this is proportional to V$(xo); the term proportional to 6x is 
zero. By an appropriate choice of coordinate system the remaining term in Eq. (|8.2|) can be 
diagonalized so that 

1 



2 



(8.3) 



For a displacement 6x = {61, 61, 61), A$ would then be given by 

A<l> ~ i(5/)2v^$(fo). (8.4) 

If ^{xq) were the electromagnetic potential then V^$(a;o) = if there are no sources at 
Xq. Since the condition for a stable equilibrium is that A$ > for any displacement 
6x, it follows that a collection of charges cannot be in stable equilibrium under purely 
electrostatic forces, which is Earnshaw's theorem. For many-body neutrino-exchange forces 
V^$(a;o) is non-zero in general, but it can be either positive or negative depending on the 
product (-1)^/2 cos(l/r/j) in Eq. ([O^ ). Since A$ > cannot be ensured for an arbitrary 



distribution of neutrons, it follows from Eq. ( |8.4| ) that in general such a distribution is 
unstable if neutrino-exchange forces are the dominant or exclusive forces present. This is 
the same situation that obtains for electromagnetism, and hence a similar conclusion follows: 
In the electrostatic case a stable equilibrium requires the presence of other (non-electrostatic) 
forces. In the present circumstance neutrino-exchange forces must be similarly balanced by 
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other forces (e.g., gravitational, strong), but for this to be the case, the magnitude of the 
neutrino-exchange force must be comparable to that of the other forces. Since the energy 
density arising from the gravitational or strong forces in any subvolume is smaller than the 
mass contained in that subvolume, we can then assume that the same must hold true for 
the contribution from neutrino exchange. 

It follows from the preceding discussion that the neutrino-exchange energy W{r) inside 
a volume of radius r < Rio must be smaller than the mass M(r) inside that volume. From 
Eqs. ( ^.241) and (|7.34| ) we note that for the a = (3 = term in Eq. (|7.35| ), the above condition 
leads to 



1 



M(r) 



4 



Gp\an\N{r)e 



N{ry 



r N(r) 



< 1, 



where cos{1/Tr) has been approximated by unity. It is convenient to rewrite Eq. 
introducing the length scale L defined by 

1 



L 



/ V2GF\an\p \ ^ 



1.7 X 10"^ cm 



1.1 eV 



15) 
by 

(8.6) 



In the approximation of neglecting the binding energy of the neutron star, so that M(r) 
N{r)mn, Eq. (|8.5|) can be rewritten in the form 



iv(r) rN'^(r) 
< 



(8.7) 



where in = h/rrinC = 2.1 x 10~^^ cm is the Compton wavelength of the neutron. Taking the 
logarithm of both sides of Eq. ( p.7|) then leads to the condition 

1 



m > 



In 



1 , (rN'^irX 
In 



.8) 



.LJ N{r) V 44 , 
Since N{r) = Anr^p/S, the right hand side can be expressed directly in terms of r/L = x, 
so that Eq. ( p.8| ) becomes 



Lm > 



X 



Ina;- — (2x 10"^^ + 5 x 10"^Mna; 
x-^ 



(8.9) 



This equation must hold for all values of r < i?io, which is the subvolume radius, and hence 
for all X < Rio/L. From Eq. ( ^.7|) it follows that when x = r / L < 1 the inequality holds 
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even when m = 0, and hence x < 1 is uninteresting. For x > 1 the coefficient of in 
Eq. ( ^.91 ) is small compared to In a; and hence the inequality becomes 

Lm>-hix. (8.10) 

X 

Since this inequality must hold for all x, Lm must exceed the largest value that (l/x)lnx 
can assume, which is 1/e. This gives 

Lm>l/e, (8.11) 

and, 

mc^> ^^^'""'^ = 0.4eV. (8.12) 

We note that m is proportional to the product Gpp which is the only relevant quantity 
having the dimensions of mass that can be formed from the available dynamical variables. 
Since the product Gpp also arises in the Mikheyev-Smirnov-Wolfenstein (MSW) mechanism 
for neutrino oscillations in matter, a few comments are in order relating the present 
work and the MSW effect. The effective energy -Ee// for a real neutrino of mass m and 
momentum p propagating through a neutron star is given by 

Eeff ^p + m^/2p + V2GFp, (8.13) 

and hence real neutrinos can be viewed for some purposes as if they had a small mass 
This heuristic picture helps to explain why the index of refraction for neutrinos differs from 
unity, in analogy to the index of refraction for light propagating in a medium. However, if we 
pursue the electromagnetic analogy we note that even in a dielectric medium electrostatic 
forces arising from the exchange of virtual photons still obey Coulomb's law, albeit with an 
attenuated strength. In the present case, the fact that neutrino exchange would retain its 
long-range character is significant since this implies that the combinatoric factor entering in 
Eq. ( |5.1| ) would remain as (^^^ . This in turn implies that unless the effective Fermi constant 
in the medium differed from the vacuum value by a factor O{10^^), the neutrino-exchange 
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energy density problem would still exist. In fact one would not expect to be significantly 
modified by the presence of a medium because there is no analog for neutrino exchange of 
a polarization charge in electromagnetism. Stated another way, there is no mechanism for 



shielding the neutrino-exchange force . 

Returning to Eq. ( ^.12| ) we note that the lower bound applies separately to i^e, ^'a*, and z/^, 
and is compatible with the upper bounds quoted in Eq. (|1.1|) for the three neutrino species, 
as shown in Fig. For v^. the upper and lower bounds are sufficiently close to suggest that 
direct evidence for vriy^ ^ could be forthcoming in the foreseeable future. Indeed it may 
well be the case that the anomalies in the flux of solar neutrinos discussed in Sec. I could be 
a signal for a non-zero neutrino mass compatible with the bound in Eq. ( ^.12| ). In addition, 
the implication that both and Vr must also be massive may help to solve the "missing 
mass" problem discussed earlier. 
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APPENDIX A: NOTATION, METRIC CONVENTIONS, AND STANDARD 

MODEL COUPLINGS 



In this Appendix we summarize our metric conventions and those for the Dirac equation. 
We have employed the Pauh metric conventions of Akhiezer and Berestetskii , deWit and 
Smith [BO], Lurie pT|, and Sakurai [B^. Reference [BO] contains an excellent summary of 



these conventions, along with tables relating the Pauli metric conventions to those of Bjorken 



and Drell |63| who use real 4-vector notation. In the Pauli conventions, the Dirac equation 



in configuration space for a particle of mass m is given hj {c = h = 1), 

(7 ■ a + m)^(x) = 0, (Al) 
where the Dirac matrices 7^ = 7^ satisfy 

{7/.,7j = 25/.^, {75,7^.1 = 0, 

75 = 75 = 71727374 = (l/4!)e^^Ap7M7i.7A7p- 
Here e^j^Ap is the completely antisymmetric permutation symbol, and the dagger denotes the 
Hermitian adjoint. In discussing the many-body contributions from systems containing an 
odd number of particles, the following trace identities involving e^^^Ap are useful (tr denotes 
the trace over Dirac indices): 

tr(7p7^7A7p75) = 4e^^Ap, 
tr(7p7<77r7p7;^7A75) = 4((5p^e^^^A - V^<^m!^a + ^arCppi^A (A3) 



where p, a, r, fx,u, \ = 1, 2, 3, 4. Other useful trace identities are given in Ref. [§0 . 



For purposes of deriving the Schwinger formula for the weak energy W in 

Appendix B, the effective low-energy Lagrangian describing the coupling of neutrinos to 
quarks and leptons is required. Using Ref. p the neutrino-quark interaction is given by 

G 

^7 = :^^/^(^) [^Lij)iqj (xh^^il + 75)qjix) + eR{j)iqj{x)-f^{l - 75)gi(a;)] , (A4) 
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^^{x) = i^p{x)j^{l + 75)^(a:), (A5) 

where iIj{x) and qj{x) are the field operators for the neutrino and for quark species j respec- 
tively. In the absence of radiative corrections the parameters eiij) and eij(j) are given in 
terms of the weak mixing angle 6^ by H 



(A6) 



^l{u) = ^ - ^ sin^ = sin^ Ow, 

e lid) = — + ^8110^ 9w eR{d) = ^sin'^ 9 w 
The neutrino-electron coupling can be similarly expressed as 

^Tix) = Mxh, {g^y^ + g^^^,) e{x)] , (AT) 

where e(x) is the electron field operator, and the constants g'^ and g^j^ are given by 

<?- = l + 2sin2^^ g7 = -\- (A8) 

Note that in (|A8|) the charged-current contribution has been included along with that from 



the neutral current 64 



For present purposes, we are interested in the coupling of neutrinos to a static source 
of unpolarized neutrons, protons, or electrons. Since the only relevant contributions in this 
circumstance are proportional to qjj4qj and 6746, the effective vector charges Oj describing 
the couplings to neutrons, protons, and electrons are 

an = 2[eL{d) + eR{d)] + [eL{u) + e/j(u)] = (A9a) 
ap = 2[eL{u) + en{u)] + [eL^d) + e^(d)] = ^ - 2 sin^ 9w = 0.036, (A9b) 
ae = gv=^ + 2 sin^ 9w = 0.964 (A9c) 

From Table 26.2 of Ref. |0 we note that the agreement between theory and experiment for 
the Standard Model couplings in Eqs. ([A^)-([A9|) is typically at the level of a few percent. It 
follows from the preceding discussion that the effective low-energy neutrino-neutron coupling 
can be written in the form given in Eq. ( pT^) below. 
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APPENDIX B: THE SCHWINGER FORMULA FOR W 



As discussed in Sec. II, our derivation of the fc-body (A; = 1, . . . , A^) neutrino-exchange 
potential utilizes the formalism developed by Hartle for the 4-body case ||18| , |29|| , which follows 
in turn from the Schwinger formula |2^, Eq. (|B34|) below. To help clarify the assumptions, 



notation, and metric conventions that underlie our results, we present here a derivation of 
the Schwinger formula due to Hartle ||2^. Other useful results can be found in Ref. 



which deals with the related question of effective Lagrangians in quantum electrodynamics. 

We are interested in computing the weak-interaction energy W of a collection of 
particles (e.g. a nucleus or a neutron star) due to neutrino exchange. Following the discussion 
in the Introduction, W can be viewed as the weak-energy analog of the static Coulomb energy 
Wc of a collection of electric charges, which for a nucleus can be approximated by the 2- 
body contribution in Eq. ( |2.6| ). As we have noted previously, one of the novel features of 
neutrino-exchange is that W is dominated by many-body contributions. We have shown in 
Sec. V that if W is expressed in the form 



N 

w = j2w^'"\ (Bl) 

k=2 

where W^'^^ is the fc-body contribution, then the dominant contributions to W arise from 
terms with k c:^ N. 

Consider for the sake of concreteness an idealized non-rotating spherical neutron star 
containing A^ neutrons. To derive the Schwinger formula for W, the weak interaction energy 
can be viewed (to lowest order in the Fermi constant Gp) as the energy difference between 
a neutrino propagating in the "vacuum" |0) inside the neutron star, and one propagating in 
the usual matter-free vacuum |0). Thus, 

W = {6\H\6) - {0\Ho\0) =S-So, (B2) 

where Hq is the free Hamiltonian for the propagating neutrino, and H includes the inter- 
actions with the neutrons. If ipi^) denotes the field operator for the interacting neutrino 
then 
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H = (fxH{x) = / d^x n{x)ip{x) = i / (fx ip'' {x)dt'ip{x), 



(B3) 



where dt = d/dt. From Eqs. (p2|) and ( p3|) , we then have {ip = ip^'j^) 

S = i j d^x (0|V^(x)74atV'(a;)|0), (B4) 

and an analogous formula for £q. We can express S in terms of the neutrino propagator by 
writing Eq. (p^ ) in the form 



= 1 f d'x{{-f^)c,pdt{6\Mx')M^)\^)}^,_ 



(B5) 

where a and (3 are spinor indices and dt acts only on unprimed variables. We can assume 
without loss of generality that t' > t in which case, 



(B6) 



If we define the neutrino propagator Sf{x,x') by 



[SF{x,x')]i3a = {6\T[iJi3{x)iJa{x')]\6) , 



then S can be cast in the form 



(B7) 



£ = -i J d^x {dttT^y^Spix, x')]} 



(B8) 



where tr denotes the trace over Dirac indices. Following Schwinger we introduce the Fourier 
transform Sf{x,x' , E) defined by 



oo 2% 



Sf{x, x') = 
Combining Eqs. (p8| ) and (p9|) we have 

J-oo 2n J J 

/ dE E^4Sf{x, X , E) 

.J —oo 



m 



1 



d X < tr 



(BIO) 
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where the hmit t' ^ t has been taken following the action of dt- If the order of integration 
with respect to x and E is interchanged then 

£ = — r dE E I d^x hihiSF{x,x,E)]\ 

= zi r dE E{i:i[^,Sp{E)]}- (Bll) 

In the last step of Eq. (|B 1 1| ) we have introduced the operator Sf{E) whose matrix elements 
give the c-number function Sf{x,x' , E): 

{x\Sf{E)\x) = Sf{x,x,E). (B12) 

The limit x' ^ x followed by / d^x can then be viewed as a formal trace in configuration 
space, so that in the notation of Eq. ( |B11| ) we can write symbolically, 

Tr = tr X J d^x (B13) 

The expression for S in Eq. (pil| ) can be recast into a more useful form by expressing the 
interacting neutrino propagator Sf{E) in terms of the free propagator Sp\E). Following 
the discussion in Appendix A, we assume that the low-energy coupling of the neutrons and 
neutrinos can be expressed in terms of the effective Lagrangian density 



Cj{x) = ^a^N,{x)i,{x). (B14) 



Here N^{x) is the neutron current, ^/^(a;) is the neutrino current, 

i^ix) = #(x)7^(l + 75)^(x), (B15) 
and a„ = —1/2. From Eqs. ( [B14|) and (|A1D the complete Lagrangian density for neutrinos 



IS 

C{x) = Cq{x) + Ci{x) 

= -^{x) [7 ■ 9 + m],lj{x) + ^a„Ar^(x)^^(x), (B16) 

where m is the neutrino mass, which we will assume to be zero at this stage. The equation 
of motion for il){x) can be obtained from the Euler-Lagrange equation, 

74 



_d_ 

dx) 



dC 

d{dip/dxx) 



dC 

dip 



and is given by 



7-a-^7-iV(l + 75) 



ij{x) = 0. 



It follows from Eq. ( |B17| ) that Sp{x, x', E) is a solution of the equation 



(B17) 



(B18) 



'-f ■ 1] — J ■ N Sf{x,x' , E) = —i6 {x — x) 



(B19) 



where 



(B20) 



7 ■ ?7 = 7 ■ 9 - 74-E, 



(B21) 



and d = d/dx. If the state vectors are normalized such that 



{x\x) = — sf) 



(B22) 



then Eqs. (pT8|) and ( pT2|) lead to 



b-V-r N]SpiE) = -tI, 



(B23) 



where / is the identity operator. Multiplying both sides of Eq. (|B23| ) by 74 gives 



[E - H]SpiE) = ^74, 



(B24) 



where 



H = —id ■ d — 747 ■ N = Hq — 747 ■ A^, 



(B25) 



using a = ^747. It follows that 



SpiE) = t[E - H]-'^,, 
sP{E)=t[E-Ho]''^,. 



(B26a) 
(B26b) 
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Combining Eqs. (pil|) and ( P26a] ) we can express £ in the form 

1 



— 7 ^OO 



S = — dE E\Tt ^ ^ 

2nJ-oo I \E-H 



(B27) 

and carrying out a partial integration allows Eq. ( |B27| ) to be written as 

dE\Yi{E - H)^ . (B28) 
Since £q can be obtained from Eq. ( P28|) by substituting Hq for H, we find for W 



E = —TilElniE - H) 
27r 



W = £-£o 

—i 



-Trl^ln 

— Tr / rfEln 

2ti U-oo \E-Ho 



°° [E-H 
dEln 

oo \E — Hq 



(B29) 



Using Eqs. ( |B24| ) and ( [B26bD the remaining term in Eq. ( |B29| ) can be written as 



^ , 747- N ~ oWfTp-, 

1 + — — = 1 - ^747- N 5>'(i?)74. 



E — Hq E — Hq 

Combining Eqs. (|B29|) and ([B30|) gives 



(B30) 



W = ^Tr dE ln[l - 2747- iV 4°^(^)74] \ ■ 



(B31) 



Since the integrand represents the infinite series 



00 



ln(l-A) = -Eir 



(B32) 



k=l 



it follows that each term in the series will be of the form 



tr[747 ■ Ni-z)sP{E) ■ ■ - 74] = tr[7 ■ N{-i)S'^>{E) ■ ■ ■] 



(0), 



(B33) 



using the cyclic property of the trace. The Dirac matrices 74 can thus be dropped from 
Eq. ( [B31| ) which then leads via Eq. ( pjl9| ) to the Schwinger formula for W [0,^: 



= ^Tr I /^^ dE ln[l + ^N,^,{1 + ^,)SP{E)]] . 



(B34) 
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Following Ref. we show in Sec. Ill that the Schwinger formula leads directly to a 



finite expression for the 2-body potential V^'^^ (^12); without having to resort to the regulariza- 
tion methods employed by either FS |T5| or HS |[T^. Evidently the Schwinger formula must 
also incorporate a regularization procedure (since the full 2-body amplitude for neutrino- 
exchange is divergent), but this regularization is built in at the outset when W is expressed 
as the difference {S — Sq) in Eq. (|B2|) . 

As noted at the beginning of this Appendix, W is the analog of the Coulomb energy Wc 
in Eq. (|2.6| ) in the sense of incorporating both the integration over the charge distribution 
and the combinatorics associated with these charges. For example, when Eq. (P34|) is 



expanded to 0{Gp) as in Eq. (|3.36|) , the integrations over all space are explicitly indicated, 
and the combinatoric factors arise from counting the number of ways that particles can 
be assigned to the coordinates Xi,...,X4 in T^^xa{xi,X2,Xs,X4). In practice we will not 
calculate W directly from Eq. ( P34|) , but rather use the Schwinger formula as a generating 
functional to obtain the A;-body potentials V^^\ri2, . . . , r^i), as in Eqs. ( p.l6| ), ( |3.50| ), ( p. 531 ), 



and ( p. 551 ). These potentials will then be integrated over the spherical volume in Section IV 
to produce the U^''\ by adapting some formalism from geometric probability. The final 
expression for U is then obtained in Sec. V by supplying the combinatoric factor , which 
allows us to write 

N N 

even even 

In order to apply the Schwinger formula in Eq. (|B34| ) it is necessary to exhibit the explicit 
functional form of the neutrino propag ator SP{E). From Eq. we see that 

SP{x,x') = (0|T[^(x)^(x')]|0), (B36) 



where our conventions for Sp^ (including various factors of i) follow those of Lurie |]6T| . 
Since Sp\x,x') is translationally invariant, we can set x' = without loss of generality in 
which case 

(x, 0) = SPix) = -(7 ■ 9 - m)AP{x) ^ -7 ■ dAPix). (B37) 
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Here A^!^"* (x) is the free propagator for a (massless) scalar field, which is given in configuration 
space by 

1 



47r^(x^ + ie) 



(B38) 



Inverting Eq. ( P9|) we then have 



dt e'^^'sPix) 



[7-(9-74S] / dt e 



JEt 



■1) 



(B39) 



47r2(r2 - t2 -j- ^g) ■ 

The integral in Eq. (|B39| ) can be evaluated using contour integration, by closing the contour 
in the upper half-plane for E > 0, and in the lower half-plane for E < 0. The result is 



sP{x,E) = j-r^ 



Arc \x\ + ie 



(B40) 



where 



7 ■ 7] = 7 ■ 9 — 74-E. 



(B41) 



To describe the propagation of a neutrino from fj to fi where fi and fj are the coordinates 



of two neutrons, let x — ^ fij = r j — fj, with = \fij\. Thus the explicit expression for 
Sp\fij, E) to be used in evaluating the many-body contributions is 

J\E\(rij+it) 



E) = [^-d,,-^,E] 



I e 



477 + ie 



= 7 ■ r]{i3)AF{fij,E), 



(B42) 



where (9jj =d/dfij. In the Tr notation of Eqs. ( |B13D and ( [B34| ), the propagator Sp\fij, E) is 
to be thought of as the ij matrix element of the operator Sp\E) in Eq. ( |B34] ). We note that 
when ln[l -!-■■•] in Eq. ( |B34| ) is expanded using Eq. ( |B32| ), the integrand in Eq. ( |B34| ) will 
contain a polynomial in E in which odd powers of E can be dropped due to the symmetric 
integration limits. 
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APPENDIX C: PROBABILITY DISTRIBUTIONS FOR POINTS IN A 

SPHERICAL VOLUME 



As noted in Sec. II, the average value {g) of a function g{r) taken over a 3- dimensional 
spherical volume of radius R is 

(CI) 



(g) = / drV3{r)g{r), 
Jo 



where r = = |ri — r*2| is the separation of two points, and Vslr) is the probability 
distribution. The functional form of 'P„(r) for an n-dimensional ball has been discussed by 



a number of authors whose work is summarized in Refs. 66 and 67 . It is convenient to 



introduce the scaled dimensionless variable s = r/2R which satisfies 1 > s > 0, and to 
re-express Vs^r), g{r), and {g) in terms of V-si^s) and g{s) so that 



{g) = dsVs{s)g{s). 
Jo 



(C2) 



Vn{s), which denotes the probability that the scaled distance between two points in an 
n-dimensional ball will be in the interval (s, s + ds) is then given by |]66| , |67 



Pnis) = 2^ns^-'h_,2[{n + l)/2, 1/2]. 

Ix{p, q) is the incomplete beta function defined by 

T{p + q) 



(C3) 



I U a) = tSP±ll [ dtt^-Hl-tY-' 
^AV.q) Y(p)T(q)Jo ^ ' ' 



(C4) 



T{p)T{q) 

The results for n = 1, 2, 3 are of particular interest in physics, and the corresponding prob 
ability distributions are given by ||66|j6 



-Piis) 



2(1 - s 
16 



cos s 



12sUl 



(2 + 3). 



(C5) 
(C6) 
(C7) 



The expression for V^{s) in Eq. ( |C7| ) has been obtained independently by Overhauser |^ 
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The following properties of P„(s) will be useful in the ensuing discussion: 

CdsVn{s) = l, (C8) 

JO 

Vn{l) = 0, n > 1; P„(0) = 0, n > 2; (C9) 

P;(l) = 0, n > 2; V'M = 0, n > 3, (CIO) 

where the primes denote differentiation with respect s. Another useful result is the mean 
value 's{n) of the separation of two points in an n-dimensional ball which is given by p6| 



-sin) - ^ r(n + 2)r(n/2) 

2\n + l) r(n + 3/2)r((n+l)/2)- ' 

For present purposes we are interested in = 'P(s), which is used in Sections II and IV 
to evaluate the 2-body contributions arising from photon-exchange and neutrino-exchange 
respectively. As can be seen from Eq. (|C7|), Vil) = V{0) = 0, as required by Eq. (|C^); 
Eq. ( |C10|) can be verified by noting that 

= 12s(l-s)(4-5s-5s2). (C12) 

Since both V{s) and V'{s) vanish at the endpoints of the physical region (s = 1 and s = 0), 
it follows that V{s) is a steeply falling function of s at the wings of the distribution. The 
behavior of V{s) near s = and 1 can be verified from Fig. |^ which exhibits V{s) and 
V'{s) in the interval 1 > s > 0. Two other quantities of interest are the location sq of the 
maximum of V{s), and the mean separation s(3) of two points in a 3-dimensional sphere. 
From Eq. ( pi2|) V'{s) = gives local minima for V{s) (in the physical region 1 > s > 0) 
at s = 0, 1 and a local maximum for V{s) when (4 — 5s — 5s^) = 0. The physical root is 
then So = [—1/2 + (21/20)^^^] = 0.5247 which is close to the middle of the physical region, 
s = 0.5. To evaluate s(3) we use Eq. ( |C11| ) with n = 3 which gives s(3) = 18/35 = 0.5143. 



The result s(3) ~ 0.5 is understandable given that V{s) peaks near sq — 0.5, and the slight 
difference between s(3) and sq is a reflection of the fact that V{s) is not symmetric about 
s = 0.5. 
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For purposes of averaging the fc-body 0-derivative contributions over a spherical volume 
it is in principle necessary to know the fc-body generalization of V{s) in 3- dimensions, and 
to be able to integrate this distribution over the coordinates of k < N = (9(10^^) particles. 
At present the dependence of this function, V^^\sij) = V^'^\si2, S23, . . . , s^i), on the vari- 
ables Sij = — rj\/2R is not known. Moreover, even if V^^\sij) were known, the task of 
evaluating the A;-body generalization of the integral in Eq. ( [4.21|) would be beyond present 
computational capabilities. For these reasons we use the preceding discussion of V{s) to 
approximate l^^^-'(ri2, ^23, • • • , ^fci) in Eq. ( p. 591) by replacing r^j with its mean value. 



{nj) = 2Rs{3) ~ R. (C13) 



The integrals in Eq. ( [4.16|) can then be carried out as in Eq. ( [4.17]) and give 



Uq — 



where rjj = 2Rsij. Using the normalization condition, Eq. ( [4.15[ ), we find. 



(C15) 



° 7iR \27tV2R^J k^+^' 
and this is the expression which will be used in Sees. IV and V. The "mean value approxi- 



mation" in Eq. ( P13|) can be justified by comparing Eq. ( |C13|) to the actual value of {rij) 



for the 2-body electromagnetic, 3-body weak, and 4-body weak potentials which have been 
evaluated directly We find for these cases, respectively: (rj^) = 0.83-R, (rjj) = 0.48i?, 



('"ij) = 0.62i?. We note that the mean value approximation in Eq. ( pi3| ) overestimates 
(rjj), and hence it underestimates |f/o^''| and ultimately \W\ in Eq. (|5.37|) . Moreover, the 
uncertainties arising from the mean value approximation are no worse than those inherent 
in estimating R itself [^,^. The fact that all the values of (r-y) are reasonably close to R 
can be understood intuitively as follows: The constraints in Eqs. ( P9| ) and ( |C1CI| ) serve to 
suppress 'P(s) at the wings of the distribution and, since V{s) is normalized, the result is a 
peaking of 'P(s) near s ^ 0.5 which corresponds to rij ~ R. 
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APPENDIX D: EVALUATION OF J2k ^^''^ 



From Eq. (|5.3|) the weak energy W is given by 



even 



In Sec. V we approximated Y,k by replacing k^^'^ in the denominator of Eq. (pi|) by N'^^"^, 
noting that the sum is dominated by terms with k c::^ N. In this Appendix we show how 
can be evaluated with the factor of k^^"^ included, should a more refined approximation to 



the sum be desired [701. Define 



and 



k k 



In Eq. (|D|) let Tr = so that 



The last step in Eq. ( p4|) can be verified by substituting t = (y — z) which allows the integral 
in Eq. (|D4| ) to be expressed in terms of the gamma function: 



, oo 

dz e^'iy - zy/^ = e'^y / dt e-^H^''' = -^r(3/2). (D5) 
Returning to Eq. ( p4[ ) we substitute z = In to recast J2k i^^^o the form 

Hence, 

Using Eq. ( |D2| ) we see that the expression for Yl,k the right-hand side of Eq. ( p7| ) now has 
the same form as the sum previously evaluated in Sec. V starting from Eq. ( |5.7| ). Hence by 
utilizing the result in Eq. ( p. 201 ), the sum on the left-hand side of Eq. ( |D7| ) can be expressed 
in terms of a one-dimensional integral, which can be evaluated numerically if necessary. 
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APPENDIX E: THE SCHWINGER-HARTLE FORMALISM FOR MASSIVE 

NEUTRINOS 



In this Appendix we generalize the Schwinger-Hartle formahsm presented in Appendix B 
to the case of massive neutrinos. Although the Schwinger formula itself retains the same form 
as in Eq. ( [B34| ), the expression for the free neutrino propagator Sp\E) is modified when 
the neutrino mass m is different from zero. The final expression for the massive propagator 
Spl^{E) is given in Eq. (|E18|) , and the application of the Schwinger-Hartle formalism to the 
present problem is discussed in Sec. VII. 

To establish our conventions when m 7^ 0, we begin with the expression for the propagator 
of a massive scalar field, which is given in our metric by ITTIJT^ 



1 r r dkn pi{k-x-koxo) 



A^Zix) = / d-'k — — -. (El) 

Here uj^ = + m^)^/^ and Cp is the usual Feynman contour in the complex k^ plane. 
The integral over k^ is straightforward and gives 

Since A^^(x) is an invariant function of "lie 3j ' 

it will have the same form for spacelike 
or timelike x^. In the former case A^^(x) can be easily evaluated by setting Xq = which 
gives a72[, 

A?i(^) = ^ / dy sinhye^'"-^'^^^ (E3) 

where r = |x| and sinhy = \k\/m. A^^ can be expressed in terms of the derivative of the 
Hankel function Hq^^ which has the integral representation ||72 | 



, , o roc 

H^^\r) = - dye-''''°'^y. (E4) 

TT J-00 



It follows from Eq. (p|) that 



H^^\-ir) = - / d^/e'"^^°'^^ (E5) 
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and hence fff2 



Syrr dr 



Sur —imr 



The expression for Aj^^(x) can be written in manifestly covariant form by letting r 
so that finally fT^ , 



57r 



It is useful to check Eq. ( ^Tj) in the two limiting cases of interest to us: a) For m 
use the series expansion ||7^ 



(E6) 



(E7) 



we 



(E8) 



to obtain 



rn^O 



m 



57rvx 



1 



(E9) 



in agreement with Eq. ( [B38| ). b) When \z\ is large HI)^\z) can be represented by the 
asymptotic expansion [|73| 

2 



exp{-z[2 - (z/ + l/2)7r/2]} 



(ElO) 



Hence with z = —imy/ we have approximately, 

.2 / 1 \3/2 



4^2 (" 



1 



(Ell) 



which shows the characteristic exponential decrease of Apl^{x). The analogous behavior for 
neutrinos eventually leads to the "saturation" of the many-body contribution to W, and 
thus to a resolution of the neutrino-exchange energy-density catastrophe. 

The massive fermion propagator in configuration space can be obtained from Eq. ( [E7|) 
and, as in the massless case, we are interested in calculating 



/OO /* c 

dte^^'sP^{x) = {-^-r^ + m) / 
-OO J— ( 



(E12) 
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where 7 ■ r/ = 7 ■ 9 — 74-E', and t = xq. We define 



2 dt cos{Et) 



im 



dt e 



2\ Hf\-imJr'^ -t^) 



iEt 



',^2X Hf\-imJr^-t'^) 



in 



57r 



(E13) 



where the last step follows by noting that the expression in square brackets is an even 
function of t. The integral in Eq. ( [E13| ) can be evaluated by making use of the relations 



72,75 



dt H'i^ {a^Jp^ - t^) cos(7t) 



i exp(— i/^vo^ + T^ 



(E14) 



dz 



(E15) 



Differentiating both sides of Eq. ( Pi)14| ) with respect to /?, and using Eq. ( Pi)15| ) we find 



Hi^\aVf3'^ -t'^) cos(7t) 



dt 



If we identify a = —im, P = —r, and 'y = E, then 

■ . _2 



exp{—i(3J a'^ + 7^). 



and 



57r / \m^r 



exp[ir\/ E"^ — m?), 



(E16) 



(E17) 



Sp ix, i?) = (7 ■ (9 — 74-E — m)- — expf^rV E"^ — m^) 
= {'y ■ r] - m)AFm{r,E). 



(E18) 



In the m — > limit the expression in Eq. ( [El 81 ) reduces to the massless result given in 
Eq. ( PI20 . 

When the Schwinger formula in Eq. ( [B42| ) is expanded in a perturbation series, and the 
Dirac traces are evaluated, integrals of the form 



/oo 
dEE"" exp{izVW 
-00 



(E19) 
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arise which are the m 7^ analogs of /^(-z) in Eq. ( |3.11|) . Since Fn{z) = when n is odd we 
can write 

2 fn~ dE exp(izy/E^ - m'^) = F„ (z) even n 
= { ' ^ ' (E20) 

odd n. 



The integrals in Eqs. ( [El 91 ) and ( |E20D can be evaluated recursively beginning with Fq{z) 



which may be cast in the form 



F^{z) = 2 dpe^^^^^=^, (E21) 



where = E'^ — ra^. Fq{z) can be expressed in terms of a modified Bessel function by 
analytically continuing the integral ITH] 



dp ^ =uKi{tu). (E22) 

VP — M 

Here Kn{z) is the modified Bessel function defined by 

K,{x) = ^t''^'Hi'\zx), (E23) 

where Hj)^^ is the Hankel function of the first kind. Substituting u = im and t = iz then 
leads to 

Fo{z) = 2imKi{mz). (E24) 
For small x the leading term in the series expansion of K^{x) is 

K^{x) - T-^{y - \)\x-\ (E25) 

and hence for small z 

1 1% 

Fo(z) ~ 2im = — , (E26) 

mz z 

in agreement with Eq. ( |3.12| ). However, the regime of interest when m 7^ is when x = mz 
is large, in which case K^{x) can be approximated by the asymptotic series 
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KJx) 



l!8x 



2!(8x)^ 



(E27) 



Combining Eqs. ([E24|) and (|E27|) then leads to 



Fn(z) 



2im 



TC 



2mz 



1 + 



+ 



(E28) 



which exhibits the characteristic exponentially falling behavior that leads to "saturation" of 
the neutrino-exchange forces. 

To evaluate Fn{z) for n 7^ we differentiate Fo{z) twice with respect to z: 
(PF,{z) 



dz^ 



F(J2) ^ r _ ^2-^ exp{izVE 

Jo 



and hence 



fF2{z)+m^Fo{z) 



F2{z) = {~ifFl^'\z)+m'Fo{z). 



(E29) 



(E30) 



Proceeding in a similar fashion, we find upon differentiating both sides of Eq. (|E29|) twice 
with respect to 2;, 

F^{z) = F^''\z) + 2m^F2{z) - m^Fo{z). 



(E31) 



Combining Eqs. ([E30|) and ( |E31|) then leads to 

F^{z) = Fit\z) - 2m^FP{z) + m^Fo(z) 



(E32) 



It follows from Eqs.( |E30| ) and (|E32|) that Fn{z) can be evaluated recursively for any even 



n 



in terms of Fq{z) and its derivatives Fq'\z), Fq^'{z),. . . , Fq''{z). 

The expressions for F2{z), -^4(2;),. . . Fn{z), can be further simplified by explicitly evalu- 
ating the various derivatives Fq^\z) that arise in the expression for Fn{z). From Eq. ( [E24D , 

F^^ (z) = 2?m"+^ir;"^ (mz) , (E33) 

where the superscript (n) denotes the nth derivative on both sides of Eq. ( [E33| ). The nth 
derivative of the modified Bessel function Ki{x) can be expressed in terms of Ko{x) and 
Ki{x) by utilizing the relations [|77| 



^7 



-Kiiz) 



(E34) 



k\'\z) = ~K,{z)--K,{z). 



(E35) 



Differentiating Eq. ( [E35| ) and using Eq. (|E34|) leads to 



Kf\z) = -K^{z) + \l + ^\K,{z). 



(E36) 



z \ z^ ^ 

Proceeding in tliis way we can express tlie ntli derivative i^'j"'' [z) in terms of Kq{z) and Ki{z). 
For purposes of calculating the 4-body contribution in Sec. VII the explicit expressions for 
Kf\z) and {z) are needed, and these are given by 



K^^\z) = -(i + ^)k,{z)-(^-+^-\k,{z) 



z z-' 



(E37) 



Kt\z) 



2 12 



z z-" 



7 4! 



+ - Ko{z)+ 1 + - + - K,{z). 



z^ z^ 



(E38) 



We note from Eq. ( |7.24| ) that when the trace of the 4-body matrix element is calculated the 
result will depend on the functions Fq'^\z), ■m'^FQ'^\z), and m'^Fo{z) which are given by 



m'^F^'^\z) = 2i 



F!i'Hz) = 2t 



4 

m 



2! 



Ko{mz) + m + m — Ki{mz 



z 

(2m^ 12m2\ / ^ 7m^ 4! \ ^ , . 



(E39a) 
(E39b) 

(E39c) 



In the m = limit, the only term which survives is the contribution proportional to m/z'^, 
and since this term is also the source of the bound on m in Eq. ( ^.121 ) we examine it in more 
detail. 

As m — i> we have from Eq. ([E25|) 



4! . X 4! 1 4! 
2im—Ki[mz) ~ 2im— = 2i—, 



z^ mz z^ ' 



(E40) 



in agreement with Eq. ( 13.14] ). As before, we are interested in the behavior of this term when 
m 7^ and z is large. Using Eq. ([E27|) we find 



Since the expression in square brackets is slowly varying (in z) compared to the remaining 
2;-dependent factors, the net effect of having m 7^ is to replace by exp{—'mz)/z^ in 
Fq^\ This result is the basis for the bound on m derived in Sec. VIII. 
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TABLES 

TABLE L Summary of the necessary conditions for the existence of a large many-body effect. 
For each interaction, "Yes" indicates that the condition is met, and "No" indicates that it is not. 
Among known interactions only neutrino-exchange meets all three conditions. See text for further 
details. 

Condition Strong Electomagnetism Weak ^ Gravity Neutrino Exchange 

Long-Range No Yes No Yes Yes 

Bulk Matter "Charge" Yes No Yes Yes Yes 

Large Coupling Strength Yes Yes Yes No Yes 

^Z'^-exchange. 



TABLE IL Parameters for two typical white dwarfs. 



Sirius B 



40 Eri B 



M ^ 
R ^ 

Ne 

W/Muc^ ^ 



l.O53(28)M0 
O.OO74(6)i?0 
6.3 X 10^*^ 
4.1 X 10^ 

-[^q(3x105'^+28-13-. 



O.48(2)M0 
0.0124(5)i?Q 
2.9 X 10^*^ 
6.6 X 10^ 

-[^q(1x10"+28-13- 



^From Ref. |51| 
^Defined by Eq. 

^Mu is the mass of the universe given by Eq. ( ^.41 ) 
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FIGURES 

FIG. 1. The 2-body neutron-neutron potential arising from neutrino-exchange. The sohd 

(dashed) lines denote neutrons (neutrinos). 

FIG. 2. Contributions to the 4-body potential energy arising from neutrino exchange. As 
before, solid (dashed) lines denote neutrons (neutrinos). Each of the diagrams (a), (b), and (c) 
is topologically different from the others, as can be seen by redrawing the graphs as shown. For 
each of these diagrams there is another that is obtained by reversing the sense of the neutrino loop 
momentum, as we show explicitly in Fig. 4. 

FIG. 3. a) Plot of the function r{s) = Vais) in Eqs. (2.12) and (C7). b) Plot of P'{s) in 
Eq. (C12). 

FIG. 4. Combinatorics for the 4-body diagrams. As before, the solid (dashed) lines denote 
neutrons (neutrinos). For k = A there are 3! = 6 topologically distinct diagrams that can be 
drawn, corresponding to the 6 possible permutations of the integers 1, 2, 3, and 4 arranged at the 
vertices of the neutrino loop. However, diagrams (a'), (b'), and (c') are obtained from diagrams 
(a), (b), and (c) respectively by reversing the sense of the neutrino loop momentum. Hence there 
are {k - l)!/2 pairs of diagrams, [(a) + (a')], [(b) + (b')], [(c) + (c')], etc. 

FIG. 5. The 3-body contributions arising from neutrino exchange. The solid (dashed) lines 
represent neutrons (neutrinos). As noted in the text, both diagrams must be included to reproduce 
Furry's theorem for the vector contribution. 

FIG. 6. Constraints on the masses of Ug, v^, and Ut. For each neutrino (or antineutrino) the 
shaded regions are excluded either by the lower bound in Eq. (8.12) or by the upper bounds in 
Eq. (1.1). 
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